053464 


DNA  4061F 


/A 


^ TIME  DEPENDENT  AIR  TRANSPORT 

O OF  RADIATION  (TDATR) 

•sC 


Science  Applications,  Inc. 

P.O.  Box  2351 

La  Jolla,  California  92038 


January  1976 


Final  Report  for  Period  8 August  1973—30  September  1975 


O 

CJ 


CONTRACT  No.  DNA  001-74-C-0027 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODES  B325075464  N99QAXPB04512  AND 
B324074464  N99QAXPB04509  H2590D. 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  D.  C.  20305 


D D C 

nsipnaaii 

MAY  3 19TB 


UNCLASSIFIED 


security  cl  ASSlF«CATlON  OF  THIS  PAGE  <»'hen  Dflla  Entered) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT’S  CATALOG  NUMBER 


PERFORMING  organization  NAME  AND  ADPRESS 

Science  Applications,  Inc.' 

P.  0.  Box  2351 

La  Jolla,  California  92038 


11  CONTROLLING  OFFICE  NAME  ANO  ADDRESS  f i ! 

Director  Jl  Ja]U»B(J76  y 

Defense  Nuclear  Agency  i ■ u.  wuwurn  m i 

Washington,  D.C.  20305  106 


1 * MON  ITORinG  AGENCY  NAME  A ADORESSf*/  dUferent  trnm  ControUtng  Office)  15.  SECURITY  CLASS  (of  this  repor^ 

UNCLASS  v 

C/O 

15«  A«;«;<Fir Arywii  nnwrCTr.RAQTMG 

• SCHEDULE  f ~r--- 


16  distribution  STATEMENT  (of  thii  Report) 

Approved  for  public  release;  distribution  unlimited. 


IT  DISTRIBUTION  STATEMENT  (of  the  ahsfmcf  entered  in  Block  20.  if  different  from  Report) 


18  SUPPLEMENTARY  NOTES 

This  work  sponsored  by  the  Defense  Nuclear  Agency  under 
RDT&E  RMSS  Codes  B325075464  N99QAXPB04512  and  B324074464 
N99QAXPB04509  H2590D. 


19  KEY  WORDS  '’Continue  <in  reverse  side  if  necesssry  and  identify  by  block  number) 

Radiation  Transport  in  Air 
TDATR  Computer  Code 
Time-Dependent  Prompt  Gamma  Rays 
^^Time -Dependent  Secondary  Gamma  Rays 


0STRACT  fCnntinue  on  reverse  side  if  necessarv  and  identify  hy  block  number) 

This  report  describes  the  parametrization  of  time-dependent  air 
transport  results  for  prompt  and  secondary  gamma  rays  resulting  in 
a computer  code  which  is  called  Time  Dependent  Air  Transport  of 
Radiation  (TDATR) . The  report  describes  the  generation  of  the 
data  base  for  prompt  gamma-ray  and  neutron- induced  secondary  gamma 
radiation,  the  parametrization  methods  and  the  command  structure 

of  the  resulting  code  which  serves  as  the  user's  guide  for  the 

code . 


DD  , 1473  EDITION  OF  1 NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


security  CLASSIFICATION  OF  TmiS  PAGE  (Ithen  Harm  Enf^W 

3 ^ 


UNCLASSIFIED 

SECUWTV  CL>.SSlF^C^T>OH  OF  THIS  P»GEflW..n  D«l«  Bnimrmd) 

. ABSTRACT  (Continued) 

Previous  efforts  of  radiation  transport  parametrization 
involved  several  versions  of  time  independent  results  under  the 
name:  ATR  (Air  Transport  of  Radiation). 


UNCLASSIFIED 


security  classification  of  this  f»AGCriFh«ft  D«(« 


SUI-IMARY 


This  report  describes  the  parametrization  of  time  depend- 
ent air  transport  results  for  prompt  and  secondary  gamma  rays 
resulting  in  a computer  code  which  is  called  Time  Dependent 
Air  Transport  of  Radiation  (TDATR) . The  report  describes  the 
generation  of  the  data  base  for  prompt  gamma  ray,  neutron  and 
neutron- induced  secondary  gamma  ray  radiation,  the  parametri- 
zation methods  for  prompt  and  secondary  gamma  rays,  as  well  as 
the  command  structure  of  the  resulting  code  which  serves  as 
the  user's  guide  for  the  code. 

Previous  efforts  of  radiation  transport  parametrization 
involved  several  versions  of  a time  independent  computer  code 
called  ATR  (Air  Transport  of  Radiation) . 


i 

1 


1 


CONTENTS 


Page 

1.  INTRODUCTION  7 

2.  GENERATION  OF  THE  DATA  BASE  8 

2.1  INTRODUCTION  8 

2.2  DATA  BASE  FOR  PROMPT  PHOTON  SOURCES  10 

2.3  DATA  BASE  FOR  SECONDARY  PHOTONS  17 

3.  PARAMETRIZATION  OF  THE  DATA  BASE  27 

3.1  PARAMETRIZATION  OF  PROMPT  GAMMA  RAYS  28 

3.2  PARAMETRIZATION  OF  SECONDARY  GAMMA  RAYS  30 

4.  TDATR  COMMAND  STRUCTURE  35 

4.1  SOURCE  SPECIFICATION  39 

4.2  GEOMETRY  SPECIFICATION  46 

4.3  OUTPUT  RELATED  SPECIFICATIONS  51 

5.  REFERENCES  56 

APPENDIX  A:  TABLES  OF  PERTINENT  PARAMETERS  57 

APPENDIX  B:  DENSITY  SCALING  OF  THE  TIME  DEPENDENT 
BOLTZMANN  TRANSPORT  EQUATION  FOR 
NEUTRAL  PARTICLES  63 

APPENDIX  C;  NUMERICAL  EVALUATION  OF  THE  CONVOLUTION 

INTEGRAL  65 

APPENDIX  D:  SAMPLE  PROBLEMS  71 


2 


LIST  OF  FIGURES 


Figure  Page 


1 Silicon  dose  rate  at  555  m and  1.11  km  of 
infinite  homogeneous  air  at  density  1.11 
mg/cm3  for  an  8-  to  10-MeV  instantaneous 

photon  source  12 

2 Local  time  of  arrival  at  detector  = (X-Z)/C  ....  13 

3 Arrival  time  distribution  at  555  m for  photons 

single  scattered  at  cosine  0.9679  15 

4 Total  flux  at  444  meters  in  air  at  density 

1.11  mg/cc  from  a 10-°  MeV  source  band  18 

5 Tissue  dose  rate  from  a prompt  gamma  ray  fission 

source  19 


6 


7 


8 


9 


10 


Neutron  infinite  air  silicon  dose  distribution 
due  to  a 14  MeV  neutron  source  with  source- 
target  separation  of  1 km  at  a density  of  1.11 
mg/ cm3  


Secondary  gamma  ray  infinite  air  silicon  dose 
distribution  due  to  a 14  MeV  neutron  source 
with  source- target  separation  of  1 km  at  a 
density  of  1.11  mg/cm^  


Neutron  infinite  air  silicon  dose  distribution 
due  to  a thermonuclear  neutron  source  with 
source-target  separation  of  1 km  at  a density 
of  1.11  mg/cm3  


Secondary  gamma  ray  infinite  air  silicon  dose 
distribution  due  to  a thermonuclear  neutron 
source  with  source- target  separation  of  1 km 
at  a density  of  1.11  mg/cm3  


Neutron  infinite  air  silicon  dose  distribution 
due  to  a fission  neutron  source  with  source- 
target  separation  of  1 km  at  a density  of  1.11 
mg/ cm3  


21 


22 


23 


24 


25. 


3 


Figure 

11 

12 

13 

14 


LIST  OF  FIGURES  (continued) 


Secondary  gamma  ray  infinite  air  silicon  dose 
distribution  due  to  a fission  neutron  source 
with  a source-target  separation  of  1 km  at  a 
density  of  1.11  mg/cm3  

Normalized  ratio  function  for  secondary  gamma 
rays  for  the  fission  source  

Time  normalization  function  for  secondary  gamma 
rays  for  the  14  MeV  neutron  source  

TDATR  problem  geometry  


Page 

26 

32 

34 

47 


4 


j 

I 

! 


LIST  OF  TABLES 


Table  Page 

1 Maximum  and  Average  Arrival  Times  for  Single 

Scattered  Photons  for  Given  Scattering  Cosine  ...  16 

2 List  of  TDATR  Commands  36 

3 Synopsis  of  TDATR  Commands  37 

A-1  Source  and  Detector  Energy  Boundaries  for  Prompt 
Gamma  Rays,  and  Detector  Energy  Boundaries  for 
Secondary  Gamma  Rays  (MeV)  57 

A-2  Souice  Spectra  for  Neutrons  58 

A- 3 Dose  Response  Functions  for  Gamma  Rays  59 

A-4  Detector  Local  Time  Boundary  Values  for  Prompt 

Gamma  Rays  60 

A- 5 Detector  Absolute  Time  Boundary  Values  for 

Secondary  Gamma  Rays  61 


1 


1 , INTRODUCTION 


This  report  deals  with  a time  dependent  version  of  the 
Air  Transport  of  Radiation  (ATR)  code  referred  to  as  TDATR. 

This  version  contains  the  results  of  the  parametrization  of 
dose  responses  for  prompt  gamma  ray  and  neutron- induced 

secondary  gamma  ray  radiation  in  infinite  homogeneous  air. 

("12  3^ 

Previous  versions  of  ATR^  ’ ’ ' contained  time  independent 
parametrizations  for  neutrons,  prompt  and  secondary  gamma 
rays,  X-rays  and  fission  product  radiation. 

The  prompt  gamma  ray  time  dependent  data  base  was  gen- 
erated using  the  MORSE  Monte  Carlo  Code,^^^  and  the  time  de- 
pendent secondary  gamma  ray  data  were  generated  using  the 
TDA^^^  time  dependent  discrete  ordinates  transport  code.  The 
data  base  is  described  in  Section  2 of  this  report. 

Section  3 describes  the  parametrization  of  six  dose 
responses  for  each  monoenergetic  prompt  gamma  ray  source  as 
well  as  six  dose  responses  for  secondary  gamma  rays  due  to 
three  neutron  source  spectra:  fission,  thermonuclear  and 
14  MeV. 

Section  4 of  this  report  is  intended  to  serve  as  a user's 
guide  to  TDATR.  It  describes  the  individual  control  commands 
that  are  needed  to  execute  TDATR.  The  commands  are  user 
oriented,  thus,  preserving  the  philosophy  of  earlier  time 
independent  versions  of  ATR. 

The  appendices  contain  pertinent  source  spectra,  dose 
responses,  time  parameters,  density  scaling  considerations, 
the  method  employed  in  solving  the  convolution  integral  and 
sample  problems. 
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2.  GENERATION  OF  THE  DATA  BASE 


2 . 1 INTRODUCTION 

This  section  describes  the  methods  used,  problems  en- 
countered, and  results  obtained  in  generating  a time  dependent 
data  base  of  radiation  transport  in  infinite  homogeneous  air. 

Crucial  to  decisions  concerning  the  methods  used  to  gen- 
erate the  data  base  and  the  form  of  the  data  base  is  its 
intended  use.  A driving  consideration  of  the  TDATR  data  base 
is  the  fact  that  the  data  must  be  scaled  to  different  densities. 

A detailed  description  of  density  scaling  of  time  dependent 
radiation  transport  is  provided  in  Appendix  B;  here,  we 
summarize  by  indicating  the  scaling  transformation.  The  time 
dependent  flux  ({)(r,t)  from  an  instantaneous  point  source  at  time 
t at  position  r in  an  infinite  nomogeneous  medixjm  of  density  p is 
related  to  the  data  base  flux  computed  in  density  Pq  by 

(})(r,t)  = K^(t)g(rQ,tQ) 

where 

K = p/Pq 
Irgl  = K|F1 
to  = Kt 

Note  that  time  is  density  scaled.  Thus,  ideally,  the  data 
base  should  contain  a relatively  fine  resolution  in  time  space 
.since  an  order  of  magnitude  decrease  in  density  results  in  a 
decrease  in  time  resolution  by  the  same  order  of  magnitude. 
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Two  time  dependent  radiation  transport  codes  were  con- 
sidered to  perforin  the  necessary  calculations  to  produce  the 
data  bases  --  the  (T^ime  Dependent  ^ISN)  Code  and  the 

MORSE  Monte  Carlo  Code.^^^  The  MORSE  Code  has  the  advantage 
that  results  can  be  computed  in  local  time  defined  as  the 
absolute  time  after  source  emission  minus  the  time  for  light 
to  travel  from  the  source  to  the  detector.  The  TDA  Code  has 
the  advantage  of  using  less  computer  time  for  providing 
detailed  spatial  and  energy  flux  descriptions. 

The  MORSE  Code  was  used  for  the  calculations  for  prompt  pho- 
ton sources.  The  TDA  Code  cannot  efficiently  resolve  the  arrival 
and  die  away  of  the  prompt  photons  arriving  at  a detector 
with  its  restriction  to  absolute  time.  A typical  time  reso- 
lution for  TDA  for  a detector  at  1 km  is  greater  than  10  shakes 
(1  shake  = 1 sh  = 10  seconds)  at  the  pulse  arrival  time. 

With  a local  time  calculation,  however,  such  as  a MORSE  cal- 
culation, the  clock  starts  when  the  uncollided  radiation 
first  reaches  the  detector  and  the  peak  can  be  resolved, 
typically,  in  tenths  of  shakes  --  a factor  of  100  improvement 
over  TDA. 

The  TDA  Code  was  used  for  the  calculation  of  time  depend- 
ent secondary  photon  transport.  The  time  dependence  of  secondary 
radiation  arriving  at  a detector  is  driven  by  the  time  and 
energy  dependence  of  the  neutron  transport.  Photons  from  in- 
elastic neutron  interactions  arrive  on  the  order  of  10~^  - 10"^ 

seconds  after  neutron  source  emission  and  capture  photons  arrive 
-3  -2 

at  around  10  - 10  seconds  and  make  important  contributions 

out  to  1 second.  It  is  necessary  to  transport  the  neutrons 
down  to  thermal  energies.  TDA  is  ideally  suited  for  this  cal- 
culation because  the  time  resolution  available  with  TDA  is 
sufficient  and  the  code  solves  for  the  transport  of  thermal 
neutrons  more  efficiently  than  by  Monte  Carlo  methods. 
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Both  data  bases  (prompt  photons,  secondary  photons) 
represent  Green's  Functions  for  the  time  dependence  of  the 
transport.  The  data  bases  were  generated  utilizing  a "delta" 
function  source  in  time.  Dose  was  computed  as  a function 
of  the  energy  distribution  of  the  source  and  a time  parameter 
(local  or  absolute  time) . If  we  denote  the  data  base  quan- 
tities by  G(E,t'),  the  time  dependent  dose  D(t)  from  a 
source  S(E,t)  is  given  by  an  integral  over  the  source  energy 
and  a convolution  in  time 


where 


D(t)  = 


G(E,t')S(E,T)dT 


t ' = t - T 

The  following  sections  describe  the  details  of  the  data 
base  generation  utilizing  MORSE  for  prompt  photon  sources  and 
TDA  for  secondary  photons  from  neutron  sources . 

2.2  DATA  BASE  FOR  PROMPT  PHOTON  SOURCES 

The  MORSE  Monte  Carlo  Code  was  used  to  compute  the  data 
base  of  time  dependent  air  transport  from  prompt  photon  sources. 
The  same  P5,  18  energy  group  cross  sections  that  were  used  for  the 
time  independent  calculations were  employed  in  the  time  de- 
pendent data  gase  generation.  Scalar  fluence  and  five  dose  re- 
sponses: tissue,  concrete,  air,  silicon  and  tantalum  were  calcu- 

lated. The  energy  group  structure  and  dose  responses  are  given 
in  Tab’e  A-3.  Separate  calculations  are  performed  for  sources 
uniformly  distributed  in  each  of  the  18  energy  bands.  The  air 

_3 

was  homogeneous  with  density  1.11  x 10  gm/cc. 


The  doses  were  scored  on  spherical  surfaces  centered 
on  the  source  point.  Local  time  was  used  to  record  the 
time  dependence.  Local  time,  is  related  to  the  absolute 

time,  T^,  after  source  emission  by 

Tl  - Ta  - R/c 

where  R is  the  radius  of  the  detector  from  the  source  and  C is  the 
speed  of  light.  Thus,  the  uncollided  radiation  for  an  instanta- 
neous source  arrives  at  T,  =0.0.  The  doses  were  resolved  into 

^ -9  -4 

6 time  bins  per  decade  from  10  to  10  sec  local  time. 

With  P5  expanded  cross  sections,  the  MORSE  Code  utilizes 
three  discrete  scattering  angles  for  each  group  to  group  energy 
transfer.  When  a collision  event  occurs  during  the  random 
walk,  the  outscatter  energy  group  is  sampled,  then  one  of  the 
three  scattering  angles  for  the  specific  incoming-outgoing 
energy  group  set  is  selected.  A ray-effect  phenomenon  was 
discovered  in  the  early  time  dose  rates  due  to  the  discrete 
angle  scattering  treatment  in  MORSE. 

At  times  <100  ns,  the  dose  rate  computed  with  the  con- 
ventional MORSE  Monte  Carlo  Code  exhibited  nonphysical  varia- 
tions when  resolved  into  six  time  bins  per  decade  as  shown  in 
Fig.  1.  The  "humps"  in  the  dose  rate  were  attributed  to  the 
discrete  angle  scattering.  Single  scattered  photons  are  the 
primary  contributor  to  the  early  time  dose  and  the  correlation 
between  the  time  of  arrival  at  a detector  of  singly  scattered 
photons  and  the  discrete  scattering  angles  causes  the  nonphysical 
dose  rates. 

The  distance  traveled  and  the  time  for  a single  scattered 
photon  arriving  at  a detector  is  pictorially  presented  in  Fig.  2. 

A photon  travels  uncollided  to  its  first  collision  site  a distance 
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Figure  1.  Silicon  dose  rate  at  555  m and  1.11  km  of  infinite 
nomogeneous  air  at  density  1.11  mg/cm^  for  an 
8-  to  10-MeV  instantaneous  photon  source. 
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X from  the  source  point.  It  then  undergoes  a scattering 
through  angle  0 and  travels  distance  y to  the  detector.  The 
local  time  at  the  detector  of  the  photon  arrival  is  (y-z)/C. 
Elementary  calculus  will  show  that  the  maximum  time  of  arrival 
occurs  when  the  uncollided  distance  is  given  by 


where  R is  the  radius  of  the  detector.  The  maximum  time  of 
arrival  is 


The  arrival  times  of  all  singly  scattered  photons  which 
scatter  at  a given  cosine  were  calculated.  The  results  of 
a typical  calculation  are  shown  in  Fig.  3.  The  arrival  times 
are  weighted  by  an  inverse  exponential  of  the  total  number  of 
mean-free-paths  traveled  by  the  photon.  The  distribution  shows 
a sharp  peak  near  the  maximum  arrival  time  (the  latest  time 
a photon  scattered  at  the  given  cosine  can  arrive  at  the 
detector) . 

The  sharp  peak  near  the  maximum  arrival  time  gives  rise 
to  the  ray  effect  (if  it  were  uniform  there  would  be  no  corre- 
lation). Table  1 gives  the  maximum  and  average  arrival  times 
for  single  scattered  photons  for  the  four  smallest  scattering 
angles  from  Group  1 (10-8  MeV) . The  first  three  arrival  times 
correspond  to  the  "humps"  seen  in  the  data.  At  later  times, 
more  and  more  contributions  come  from  multiply  scattered  photons 
and  the  ray  effect  is  hidden. 
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The  MORSE  Code  was  modified  to  sample  the  scattering 
cosine  from  the  continuous  Klein-Nishina  distribution  for 
the  first  scattering  only  and  then  continue  with  the  discrete 
scattering  angle  selection  on  subsequent  collisions.  The 
results  of  this  computation  were  compared  in  Fig.  1 with  the 
previous  results.  The  first-collision  Klein-Nishina  sampling 
removed  the  correlation  which  gave  the  non-physical  ray  effect 
with  the  conventional  MORSE  calculation. 

Figure  A compares  the  flux  at  444  meters  computed  using 
the  Flair  Code^^^  with  the  ATR  data  base.  There  is  good  agree- 
ment between  the  two  sets  of  data.  Figure  5 shows  the  tissue 
dose  rate  at  several  ranges  for  a prompt  gamma  ray  fission 
source . 

2 . 3 DATA  BASE  FOR  SECONDARY  PHOTONS 

The  TDA  (Time  Dependent  ANISN)  Code  was  used  to  generate 
the  secondary  photon  data  base.  The  TDA  Code  determines  the 
solution  in  time  space  by  using  the  ANISN  procedure  to  solve 
for  the  flux  at  each  time  step  utilizing  the  flux  extrapola- 
tion from  the  previous  time  step  plus  any  source  emission 
during  the  time  interval. 

Proper  choices  of  spatial  and  time  meshes  are  crucial  to 
d good  solution  from  TDA,  since  they  cannot  be  selected  inde- 
pendently. A problem  involving  an  instantaneous  source  (delta 
function)  will  have  a "wave"  of  uncollided  and  few  within  group 
scattered  particles  which  move  through  the  system  at  a speed 
corresponding  to  the  speed  of  the  average  group  energy.  A fixed 
spatial  mesh  which  is  fine  enough  to  cover  the  wave  at  all 
points  in  the  system  would  require  a prohibitive  amount  of  core 
storage  and  result  in  very  inefficient,  time-consiiming  calcu- 
lations, to  iterate  through  the  fine  mesh  as  the  wave  front 
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10'^  10‘®  10'^  10"^  10"^ 


Local  Time  (Sec) 

Figure  A.  Total  flux  at  444  meters  in  air  at  density 
1.11  mg/cc  from  a 10-8  MeV  source  band. 
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A time  dependent  spatial  zoning  package  was  incorpor- 
ated into  the  TDA  Code  to  peirmit  the  use  of  fine  zones  near  the 
wave  front  (where  most  interactions  are  taking  place)  and  rela- 
tively coarse  zones  elsewhere.  At  each  time  step  the  system 
is  rezoned  to  accomodate  the  movement  of  the  wave.  Addition- 
ally, the  rezoning  is  performed  such  that  certain  ranges  input 
by  the  user  corresponding  to  detector  distances  are  always 
zone  centered.  Input  to  the  zoning  package  includes  the  upper 
and  lower  bounds  on  the  wave  velocity,  the  number  of  fine 
zones,  and  the  number  and  range  of  detectors. 

Three  source  distributions  were  used  for  the  TDA  calcula- 
tions; a weaponized  fission,  a thermonuclear,  and  a 14  MeV 
neutron  source.  These  source  distrubtions  are  identical  to 
the  source  options  in  the  ATR-4'  code  and  are  tabulated  in 
Appendix  A.  The  time  bins  and  detector  distances  are  also 
provided  in  Appendix  A. 

Typical  results  of  the  TDA  calculations  are  shown  in 
Figs.  6 thru  11.  Neutron  and  secondary  photon  dose  rates 
at  one  kilometer  are  shown  for  each  source  distribution.  The 
neutron  data  base  was  not  parametrized  in  TDATR,  only  some  of 
the  data  is  presented  here  for  the  sake  of  completeness. 
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Figure  6.  Neutron  infinite  air  silicon  dose  distribution 
due  to  a 14  MeV  neutron  source  with  source- 
target  separation  of  1 km  at  a density  of  1.11 
mg/ cm3 . 
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Absolute  Time  (Sec) 


Figure  7. 


Secondary  gamma  rav  infinite  air  silicon 
dose  distribution  due  to  a 14  MeV  neutron 
source  with  source- target  separation  of 
1 km  at  a density  of  1.11  mg/cm^. 


Silicon  Dose  (Rads-Cm'^/Sec) 


Absolute  Time  (Sec) 


Figure  8.  Neutron  infinite  air  silicon  dose 

distribution  due  to  a thermonuclear 
neutron  source  with  source-target 
separation  of  1 km  at  a density  of 

1. 11  mg/cm3 . 


3.  PARAMETRIZATION  OF  THE  DATA  BASE 


The  data  base  that  was  used  for  the  parametrization 

of  the  prompt  gamma  rays  consisted  of  data  at  22  ranges  out 

2 

to  about  390  gm/cm  and  25  local  time  mesh  points  out  to 
10~^  seconds  for  the  following  six  responses: 

1.  Total  Fluence 

2.  Henderson  Tissue  Dose 

3.  Concrete  Dose 

4.  Air  Dose 

5.  Silicon  Dose 

6.  Tantalum  Dose 

There  are  detector  response  values  in  each  of  the  18  mono- 
energetic  source  groups.  The  dose  factors  and  the  source 
energy  breakdown  are  detailed  in  Appendix  A. 

For  the  secondary  gamma  rays  there  were  three  standard 
neutron  sources  that  were  used  and  the  same  six  responses 
were  determined  for  each.  The  source  spectra  were: 

1.  Fission 

2.  Thermonuclear 

3.  14  MeV 

The  source  weight  values  are  detailed  in  Appendix  A for  the 
fission  and  thermonuclear  source;  the  14  MeV  source  represents 
the  data  in  the  12.2  - 15  MeV  neutron  source  band.  The  results 


2 

were  determined  for  85  ranges  out  to  about  400  gm/cm  . There 
were  40  time  mesh  points  in  absolute  time  out  to  the  maximum 
value  of  two  seconds. 

3.1  PARAMETRIZATION  OF  PROMPT  GAMMA  RAYS 

The  data  was  displayed  through  various  forms  of  plots 
both  as  a function  of  time  as  well  as  range  in  order  to  eval- 
uate trends.  Since  data  generated  by  Monte  Carlo  techniques 
result  in  statistical  variations,  it  was  difficult  to  ade- 
quately categorize  the  information  and  some  averaging  through 
curve  fitting  was  necessary.  There  were  two  categories  of 
data;  one  included  the  total  fluence,  Henderson  tissue  and 
tantalim  dose  values  and  the  other  contained  the  remaining 
three  dose  responses.  The  following  ratio  function  was  calcu- 
lated for  each  of  the  responses  in  the  first  category: 


D,-  (E.  .tu.Pp) 

Ri(Ei.tk,P^)  = d'Se’^  t,  p - 


. j j. 


= 1.9 


where  D . (E, , t,  , p .) represents  the  response  values  (total  fluence 

1.  J K i6 

or  dose)  as  a function  of  source  energy  band  E. , time  t,  in 

2 J 

seconds,  range  Pj^  gm/cm  . The  index  i has  values  of  1,  2 
or  6 representing  the  total  fluence,  Henderson  tissue  dose  or 

40 

tantalum  dose  respectively.  The  index  j ranges  from  1 to  18 
representing  one  of  the  18  monoenergetic  source  groups.  The 

index  Jq  in  the  denominator  is  one  for  j = 2,3 8 and  nine 

for  j = 10 , 11 , . . . , 18 . That  is,  16  ratio  functions  are  generated 
for  each  response  such  that  source  groups  two  through  eight 
are  divided  by  the  first  (highest  energy)  source  group  and 
source  groups  10  through  18  are  divided  by  the  ninth  source 
group.  The  indices  k and  I appear  simply  to  indicate  that  t^ 
and  are  time  and  range  parameter  variables. 
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The  ratio  function  was 
the  following  relationship: 

Hn 

where  the  a^^^  are  the  coefficients  of  the  parameterization  and 
the  index  I has  been  removed  from  the  range  variable  to  indi- 
cate that  it  is  being  treated  as  a continuous  variable.  Note 
that  indices,  i,  j and  k are  being  treated  as  parameters. 

The  denominator  of  the  ratio  function  has  a considerably 
more  complicated  form  and  it  was  parameterized  separately 
according  to  the  following  relationship: 

in  = bQ  + bj^T  + /^TJ+cp7^ 

+ b^T  ^ + b^T  ^ 

where 

T = )ln  t,  = -10.82,  C2  = 10.0,  C3  = 13.0 

the  bjjj  are  the  coefficients  of  the  parameterization,  Jq  = 1,9, 
and  the  index  k has  been  removed  from  the  time  variable  in 
order  to  indicate  that  it  is  treated  as  a continuous  variable. 

The  ratio  function  is  evaluated  for  those  sources  that 
require  it  and  multiplied  by  the  denominator  of  the  ratio 
function.  Then,  in  order  to  generate  D^(Ej,t,p)  for 

^ ^k+1  ^k  - '^k-t-1’  values:  . P^) 

^i^^j'Vl'Pk^'  ^i^^j’^k-Pk+l^-  ^i^^j-^k+l-Pk+l^  evaluated 
and  a double  interpolation  scheme  is  used. 


then  curve  fitted  according  to 


= aQ  + a^p  + a2P^ 
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For  the  remaining  three  dose  responses;  concrete,  air 
and  silicon  dose,  the  following  ratio  function  was  generated: 


r.  (E. 
1 J 


DjCE.  ,tk,P)^) 
^2  ’ ^k  ’ 


i = 3,4,5 


where  ^2^^j’^k’Pi^  represents  the  Henderson  tissue  dose  response. 
This  ratio  function  has  a rather  straightforward  behavior 
in  that  it  was  represented  by  an  averaged  constant  or  a series 
of  straight  lines  as  a function  of  the  range  parameter.  Thus, 
only  one  or  two  coefficients  are  needed  to  be  stored  in  order 
to  regenerate  the  ratio  function  which  is  then  multiplied  by 
D2  (Ej  , t^ , Pj^)  to  regenerate  the  D^(Ej  , tj^  , p which  are  needed 
for  the  evaluation  of  D^(Ej,t,p)  as  previously  described. 

3.2  PARAMETRIZATION  OF  SECONDARY  GAMMA  RAYS 


For  secondary  gamma  rays  there  was  a greater  similarity 
between  the  various  responses  than  for  prompt  gamma  rays,  thus, 
the  responses  were  more  readily  categorized. 

First,  the  total  fluence  was  parametrized  for  all  three 
sources  by  the  following  function: 


£n  ’ ^k  ’ P^J 


+ a,  p 


+ a2P' 


+ a3P 


^4 

+—  + 
P 


a^ 


where  D3(Sj,tj^,p)  represents  the  total  fluence  response  for 
source  spectra  S ^ ( j = 1,2,3  representing  the  fission,  thermo- 
nuclear and  14  MeV  neutron  sources  respectively)  , the  tj^  are 
the  time  parameters,  pis  the  continuous  range  variable  and  the 
a^^  are  the  coefficients  of  parametrization. 
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Then,  for  the  remaining  dose  responses  the  following 
ratio  was  generated: 


> . 

J 


Di(S.  ) 


DTTST7t^Tp^~T 

o j 


i=2,3, ... ,6;  j=l,2,3 


= R:(s..t^.P,)  N.(s..tj^) 


where  Si  = 1 . 
o 

Figure  12  represents  examples  of  j > , P for  the 

fission  source  as  a function  of  range  with  time  being  the 
parameter. 


The  ratio  function,  (S^  , tj^  , was  parametrized  accord- 

ing to  the  following  scheme.  For  i=2,  that  is  the  Henderson 
tissue  dose,  the  following  function  was  used: 

£n  " ^0  ^3^^ 

where  the  a^  are  the  coefficients  of  the  parametrization. 

For  the  remaining  dose  responses  the  ratio  function  was  treated 
as  a perturbation  of  R2  (Sj  , tj^ , p)  , that  is: 

Ri(S.,tk,p)  = R2(S^,tj^,p)  P(S^,tj^,p)  , i=3,4,5,6 

where  P(Sj,tj^,p)  is  a perturbation  fvmction  of  the  range 
variable  that  is  represented  by  simple  straight  line  functions 
cf  p . 
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The  second  term  of  the  ratio  function,  N.(S.,t,),  can  be 

1 J iC 

viewed  as  a normalization  function  in  time,  and  it  ensures  that 
the  value  of  the  ratio  function  is  1.0  at  the  first  range  value. 
Figure  13  shows  a representative  sample  of  the  time  ratio  function 
for  the  14  MeV  source  as  a function  of  time  for  two  dose  responses 
This  term  was  parametrized  separately  by  the  following  function 
of  time:  let  x = K,n  t 


£n|^N^(S^  ,t)j  = ag  + aj^T  + 


a2T 


+ a3T^  + a^x'^  for  t < 1.33  x lO"^  sec 


= bo  + b^x 


for  t > 1.33  X 10' 


sec 


where  the  a and  b are  the  coefficients  of  parametrization. 
m n 

Thus,  the  value  of  (S^  , tj^ , p^)  can  be  obtained  by; 

R . (S  . , t,  iPn)  (S.,  t.  , p « ) 

^ ^ 


and  a double  interpolation  scheme  is  used  to  find  the  response 
value  at  an  arbitrary  time  and  range  value. 
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4.  TDATR  COMMAND  STRUCTURE 


The  TDATR  command  structure  is  similar  to  that  of  ATR 
with  minor  modifications  and  additions.  A typical  TDATR 
command  is  of  the  following  general  form: 

*<COMMAND  WORD>,<UNITS  DEFINITION>,<LIST  OF  VALUES> 

All  commands  must  begin  with  an  asterisk.  If  a TDATR  command 
is  too  long  to  fit  on  a single  80-character  line,  the  command 
may  be  continued  on  subsequent  lines  not  beginning  with  an 
asterisk.  The  restrictions  on  continuation  are  that  no  single 
part  of  a TDATR  command,  such  as  a number,  may  itself  be  split 
into  two  lines  and  there  should  be  at  least  two  numbers  (values) 
on  the  first  card  image. 

The  command  word  is  a mnemonic  name  for  the  type  of  action 
to  be  taken.  Examples  of  command  words  include: 

N-SVAL  - Source  values  for  neutrons 

LTIME  - Specifies  local  time  for  output 

TITLE  - Title  of  a TDATR  problem  output  with  the  string 
that  follows  the  first  blank  after  the  command 
word. 

A complete  list  of  the  commands  for  TDATR  is  contained  in 
Table  2.  The  upper-case  letters  in  the  figure  indicate  specific 
syntactic  elements  while  the  lower-case  letters  indicate  para- 
metric fields.  Table  3 contains  a short  synopsis  of  the  individual 
commands . 

The  second  field  of  the  general  command,  the  units  definition, 
is  delimited  by  commas.  The  units  definition  field,  which  is 
optional,  serves  to  explain  the  meaning  of  the  numbers  in  the 
list  of  values  that  follows.  Typical  units  include: 


35 


Table  2.  List  of  TDATR  commands. 


1.  *N-SOURCE  values 

2.  *G-SOURCE(i) 

3.  *G-EVAL,  units,  values 

4.  *G-SVAL(i^  ±2  •••)•  units,  values 

5.  *G-STVAL,  units,  values 

6.  *G-GAUSS  to  B FWHM 

7.  *G-TIME,  units,  values 

8.  *y-NORM  value 

9.  *y-YIELD  value 

10.  *xx,  linits,  value (s) 

11.  *GR01JND,  units,  value 

12.  *RESP/z/(ij^  ±2  • • •) 

13.  *LTIME/z/,  units,  values 

14.  *ATIME/z/,  units,  values 

15.  *TITLE  n 

16 . *EXC 

17.  *ST0P 

18.  *FIN 


y = N or  G which  stands  for  either  neutrons 
or  gamma  rays . 

z = G or  NG  which  stands  for  either  gamma  rays 
or  secondary  gamma  rays. 
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Table  3.  Synopsis  of  TDATR  commands. 


I 

h 

I 

i 

I 


SOURCE  - indicates  that  the  ith  source  option  will  be  chosen 
for  prompt  gamma  rays  or  gives  the  weighting  values 
for  neutron  sources 

EVAL  - indicates  that  the  energy  group  boundaries  will 

follow 

SVAL  - indicates  that  the  energy  dependent  source  intensity 

will  be  entered 

STVAL  - indicates  that  the  time  dependent  source  intensity 

component  will  be  entered 

GAUSS  - parameters  that  define  a Gaussian  time  dependent 

source 

TIME  - time  bins  used  for  time  dependent  source  specification 

NORM  - source  normalization  value  follows 

YIELD  - yield  of  source 

XX  - indicates  the  geometry  component  and  associated 

values  will  follow 

GROUND  - specifies  the  ground  evaluation 

RESP  - a response  function  selection  for  secondary  or 

prompt  gamma  rays  will  follow 

LTIME  - the  local  time  bins  of  interest  for  time  dependent 

radiation  environment  will  follow 

ATIME  - the  absolute  time  bins  of  interest  for  time  dependent 

radiation  environment  will  follow 

TITLE  - alphanumeric  character  string  to  be  used  as  a title 

will  follow 

EXC  - an  action  command  to  execute  the  problem  specified 

by  prior  input 

STOP  - indicates  the  end  of  a problem 

FIN  - indicates  the  end  of  a TDATR  session 


MEV 
KEV 
PER  SEC 
NANO 


energy  values 
energy  values 

source  normalization  values 
time  values 


Specific  sets  of  unit  definitions  are  appropriate  for 
each  command,  as  will  be  indicated  in  subsequent  sections  of 
this  report.  A default  unit  definition  has  been  selected  from 
the  set  appropriate  to  each  command,  and  is  used  whenever  the 
units  definition  with  its  surrounding  commas  is  omitted. 

The  list  of  values  element  of  an  TDATR  command  is  used  to 
specify  the  numerical  data  that  a command  may  require.  A 
number  in  a list  may  appear  in  a variety  of  forms  to  suit  the 
particular  user  or  problem.  For  example,  some  of  the  forms 
in  which  the  number  400.  may  appear  are  as  follows: 


400  400.  4.E+2  4E+2  4. +2  4+2  4000-1 


At  least  one  number  must  appear  in  a list  of  values  element. 

Two  or  more  numbers  are  separated  from  one  another  by  the 
occurrence  of  one  or  more  blank  characters.  Therefore,  the 
user  is  restricted  from  specifying  a ntimber  in  which  internal 
blanks  appear,  or  which  is  split  on  two  or  more  card  images. 

A further  restriction  exists  upon  the  magnitude  of  such  nvimbers: 
since  a number  is  interpreted  as  a function  of  up  to  three 
integer  parts  (a  whole  part,  a fractional  part,  and  an  exponent), 
none  of  the  parts  of  a number  may  exceed  in  magnitude  the 
greatest  integer  value  appropriate  to  the  host  machine,  nor  can 
the  number  generated  from  these  three  parts  exceed  the  host 
machine's  allowable  range  of  representable  numbers. 

Section  4.2.1  contains  a description  of  other  convenience 
features  associated  with  the  entry  of  a list  of  values. 
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4.1  SOURCE  SPECIFICATION 


4. 1 . 1 General  Source  Input  Description 

The  specification  of  the  source  distribution  for  TDATR 
is  complicated  by  the  fact  that  the  prompt  and  secondary  gamma 
rays  were  treated  rather  differently  during  the  parametrization 
and  the  source  specification  must  reflect  those  differences  but 
at  the  same  time  must  also  adhere  to  the  general  framework. 

Since  there  are  only  three  possible  neutron  source  selections 
for  which  secondary  gamma  ray  data  were  parametrized,  there  is  no 
energy  dependence  as  such  but  a mix  of  fission,  thermonuclear 
and  14  MeV  sources  is  possible  through  the  use  of  the  *N-SOURCE 
command  which  allows  for  the  specification  of  the  three  relative 
weights.  The  detector  response  values  corresponding  to  the 
three  source  spectra  are  evaluated  and  multiplied  by  the  cor- 
responding weights  and  added  together.  The  weights  can  be 
normalized  by  the  use  of  the  *N-NORM  command  which  specifies 
the  total  normalization  in  neutrons/kT.  The  weapon  yield  can 
be  specified  by  the  *N-YIELD  command  which  allows  the  input  of 
the  yield  in  KT  of  the  weapon.  There  is  no  time  dependence  of 
the  neutron  source,  the  weight  factors  act  as  a delta- function 
source . 

The  prompt  gamma  ray  source  specification  allows  for  a 
full  two-dimensional  source  spectrum,  i.e.:  specification  as 
' a function  of  time  as  well  as  energy  is  possible.  The  source 

component  entry  as  the  function  of  energy  is  facilitated 
by  the  *G-SVAL  command  and  its  corresponding  structure, 
if  different  from  the  data  base,  is  specified  by  the  *G-EVAL 
command.  There  are  four  different  ways  that  the  source  values 
can  be  defined  as  a function  of  time: 
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1.  Delta  function 

2.  Gaussian 

3.  Separable  time-energy  functions 

4.  Individual  source  values  entered  into  the 

time-energy  source  matrix 

The  differentiation  of  the  various  types  is  specified  by 
the  *G-SOURCE(i)  command  where  i is  the  index  of  the  source 
types  corresponding  to  one  of  the  above  choices.  If  S(E,t) 
represents  the  source  spectrum  as  the  function  of  energy  and 
time,  then  the  energy  dependent  delta- function  values  (at  t=0) 
are  specified  by  a special  form  of  the  *G-SVAL  command. 


There  is  a special  time  dependent  shape  of  the  source  spec- 
trum that  can  be  entered  via  the  use  of  the  *G-GAUSS  command  which 
is  of  the  following  form: 


f(t) 


B e 


(t-t^) 


2o' 


2 


where  t is  the  mean  of  the  Gaussian,  B is  the  amplitude  at 
o 

t=t^  given  by.  B = (a/o  /2tt  ) where  a is  a constant  and  a is  the 
variance  that  is  related  to  the  full  width  at  half  maximum  (FWHM) 
by  the  formula; 

FWHM 

o = _ - = FWHM*.  42466 

2 n /-iln  .5 

The  *G-GAUSS  command  contains  the  value  of  t in  units  of 

o 

shakes,  B in  units  of  MeV/second  and  FWHM  in  units  of  shakes. 
Then  TDATR  will  evaluate  o and  will  place  the  time  dependent 
values  into  the  array  S(E,t)  by  weighting  the  Gaussian  values 
by  the  energy  dependent  source  values  entered  via  the  *G-SVAL 
command.  Also,  the  value  of  B will  be  modified  by  the  time 
independent  energy  spectrum  from  MeV/sec  to  photons /sec  by 
using  the  average  energy.  The  time  values  at  which  the  S(E,t) 
are  evaluated  is  entered  via  the  *G-TIME  command  which  also 
applies  to  the  third  and  fourth  type  of  prompt  gamma  ray  source. 
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The  third  type  of  prompt  gamma  ray  source  assumes  that 
the  source  spectrum  is  separable  in  time  and  energy.  That  is, 
the  source  function  S(E,t)  can  be  written  as;  S(E,t)  = g(E)f(t). 
The  function  g(E)  is  strictly  a function  of  energy  and  is 
entered  via  the  *G-SVAL  command;  the  energy  structure  is  given 
by  the  *G-EVAL  command  or  it  defaults  to  the  internal  structure. 

The  function  f(t)  is  strictly  a function  of  time  and  is  entered 
via  the  *G-STVAL  command  with  values  corresponding  to  source 
time  boundary  values  entered  via  the  *G-TIME  command.  The 
S(e,t)  matrix  will  be  filled  in  by  the  point  by  point  product 
of  the  two  functions. 


The  fourth  type  of  source  definition  allows  for  an  entry 
of  every  matrix  element  of  S(E,t)  by  the  use  of  the  *G-SVAL 
command  that  allows  for  a time  index  to  be  specified  where  the 
energy  dependent  values  are  to  be  entered.  Again,  the  energy 
and  time  structure  can  be  entered  via  the  *G-EVAL  and  *G-TIME 
commands  respectively. 

For  prompt  gamma  rays  the  commands  *G-NORM  and  *G-YIELD 

* 

commands  also  apply;  the  application  of  normalization  and 
yield  will  take  on  the  following  form.  If  N is  the  value 
of  the  normalization  in  photons /KT,  Y is  the  value  of  the 
yield  in  KT  and  S'(E,t)  is  the  source  matrix  in  per  MEV  per 
second,  then  the  final  source  spectrum  is  described  by; 


S(E,t) 


NY  S' (E,t) 
ff  S' (E,t)  dE  dt 


where  the  double  integration  becomes  a double  sum  since  the 
energy  and  time  boundaries  are  discrete  values  corresponding 
to  the  discrete  values  of  S'(E,t). 

The  maximum  number  of  energy,  time  or  source  values 
allowed  is  arbitrarily  limited  to  50. 
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4.1.2  Specific  Source  Commands 
4. 1.2.1  *N- SOURCE  values 


This  command  specifies  the  three  source  weight  values 
separated  by  one  or  more  blanks  that  will  weight  the  results 
of  secondary  gamma  rays  from  the  neutron  fission,  thermonuclear 
and  14  MeV  sources  respectively.  If  less  than  three  values  are 
entered, then  the  remaining  weights  will  be  zero.  There  is  no  unit 
definition  for  these  values  since  they  are  simply  used  as  weighting 
factors  for  the  respective  source  contributions.  Aside  from  the 
optional  *N-NORM  and  *N-YIELD  commands , this  is  the  only  command 
necessary  for  the  neutron  source  definition. 

Examples : 

*N-SOURCE  .25  .25  .5 

which  specifies  a weighting  factor  of  .25  for  both  fission 
and  thermonuclear  sources  and  a weighting  of  .5  for  the 
14  MeV  source. 

*N- SOURCE  001 

which  indicates  that  only  the  14  MeV  source  is  used. 

4. 1.2. 2 *G-S0URCEa) 

This  command  specifies  the  intended  source  spectrum  struc- 
ture of  a prompt  gamma  ray  source.  The  value  of  i is  1,2,3  or  4 
denoting  one  of  the  four  types  of  sources  that  were  detailed  in 
Section  4.1.1.  Other  accompanying  commands  will  specify  the 
value  details  of  the  source  spectrum.  Note  the  similarity 
between  this  command  and  the  one  described  in  Section  4. 1.2.1. 

Example : 

*G-S0URCE(2) 

which  specifies  that  a Gaussian-type  source  will  be  used 
for  the  prompt  gamma  ray  source  spectrum  and  the  para- 
meters of  the  Gaussian  will  be  supplied  by  the  accompany- 
ing *G-GAUSS  command. 
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4. 1.2.3  *G-EVAL,  units,  values 


units  = MEV , KEV . If  ".units,"  is  not  specified  then 
the  default  is  MEV. 

values  = energy  point  or  histogram  values  from  low 
energy  to  high  energy  ascending  order. 

Note  that  if  units  are  not  specified  than  the  surround- 
ing commas  must  also  be  omitted.  The  distinction  between  point 
or  histogram  values  will  be  made  by  TDATR  by  counting  the  number 
of  values  entered  with  a corresponding  *G-SVAL  command.  If  the 
number  of  values  are  the  same  then  a point  value  is  assumed,  if 
there  is  one  more  energy  value  than  source  value  then  a histo- 
gram value  structure  is  assumed,  otherwise  the  code  will  indicate 
an  error  condition. 

If  this  command  is  not  present  then  the  code  assumes  the 
internal  18  group  structure  listed  in  Appendix  A. 

Example : 

*G-EVAL  .01  .03  .1  .4  .8  1.2  1.8  2 4 6 8 10 

which  represents  the  energy  values  for  a prompt  gamma  ray 

source  in  units  of  MeV  which  is  the  default. 

4. 1.2. 4 *G-SVAL(i^  ±2  ■■■),  units,  values 

i^ , i2  , ■ . ■ , ij^  = ordinal  indices  of  time  bins  into  which 
the  values  are  to  be  placed  separated  by  one  or  more  blanks.  This 
construct  will  only  be  present  when  the  fourth  prompt  gamma  ray 
source  specification  is  desired  because  the  delta  function, 
Gaussian  and  separable  source  functions  do  not  have  an  explicit 
time  index  associated  with  them.  If  this  construct  is  not  pre- 
sent then  the  surrounding  parentheses  must  also  not  be  present. 

units  = PER  MEV,  PER  KEV,  PER  GROUP;  the  default  is 
PER  GROUP.  Note  that  point  source  values  with  a PER  GROUP  unit 
specification  make  no  sense  and  will  result  in  an  error  con- 
dition. 
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values  = energy  dependent  source  values  that  will  be  used 
to  compute  the  actual  source  values.  These  values  correspond 
to  a given  set  of  energy  values  provided  that  the  *G-EVAL  com- 
mand is  also  present,  otherwise  the  values  correspond  to  the  inter- 
nal 18  group  energy  structure.  These  values  must  correspond  to 
the  energy  values  in  an  order  which  is  low  to  high  in  ascending 
energy  order. 

If  the  fourth  source  option  is  used  then  perhaps  several 
commands  will  be  necessary  to  fill  the  whole  source  matrix.  The 
source  values  corresponding  to  unspecified  ordinal  indices  will 
be  set  to  zero: 

Examples : 

*G-SVAL  (1  3 5).  PER  ME V,  .01  1.5  2 2.2  .9  .02 

which  has  the  effect  of  placing  the  six  low  energy  source 
values  as  a function  of  energy  into  time  bins  1,  3 and  5. 

*G-SVAL,  PER  GROUP,  .1  1 2.2  0.5  1-3  2-4 

which  specifies  the  energy  weights  for  the  six  low  energy 
bins  for  any  one  of  the  first  three  types  of  sources. 

Note  that  the  units  definition  is  PER  GROUP  which  is 
superfluous  since  it  is  the  default  definition. 

4. 1.2.5  *G-STVAL,  units,  values 

units  = one  of  PER  GROUP,  PER  NANO,  PER  SHAKE,  PER  MICRO 
PER  MILLI,  PER  SEC;  the  default  unit  is  PER  SHAKE. 

values  = time  dependent  source  values  for  the  separable 
source  spectrum  case;  i.e.  values  of  f(t)  as  described  in 
Section  4.1.1.  The  order  of  the  values  is  assumed  to  be  from 
low  to  high  corresponding  time  values  entered  via  the  *G-TIME 
command. 

Example ; 

*G-STVAL,  PER  SEC,  2 1.5  1 .3  .1  .01  .001 

which  specifies  seven  values  corresponding  to  a time  bin 
structure  that  is  entered  via  the  *G-TIME  command. 
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4.  1.2.6  *G-GAUSS  t B FWHM 
o 

The  values  of  t , B and  FWHM  which  are  the  parameters 

of  the  Gaussian  time  dependent  source  function  as  detailed  in 

Section  4.1.1.  The  units  of  t and  FWHM  are  shakes  and  B 

o 

has  energy  units  of  MeV/sec.  The  code  will  evaluate  the  average 
energy  from  the  time  independent  spectrum  and  modify  B to 
photons/ sec. 

Example : 

*G-GAUSS  100  1 10 

which  defines  a Gaussian  time  dependent  source  shape 
centered  at  100  shakes,  has  an  intensity  of  1 MeV/sec 
and  is  10  shakes  at  FWHM. 

4. 1 . 2 . 7 *G-TIME,  units,  values 

units  = one  of  NANO,  SHAKE,  MICRO,  MILLI , SEC;  the 
default  unit  is  SHAKE. 

values  = time  histogram  boundary  values  corres- 
ponding to  the  time  dependent  dimension  of  the  source  speci- 
fication. The  values  must  be  entered  from  low  to  high. 

Example : 

*G-TIME  0 10  21.5  50  100  400  800  1+3  3+3 

5+3  1+4  2+4  6+4  1+5  1+6  1+8 

which  specifies  the  time  boundary  values  for  the 
prompt  gamma  ray  source  in  default  units  (shakes). 

4. 1.2. 8 *y-N0RM  value 

y = N or  G representing  a neutron  source  or  prompt 
gamma  ray  source  normalization  respectively. 
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value  = total  source  normalization  in  particle/KT.  There 
is  no  default  value  for  total  source  normalization  and  if  this 
command  is  not  specified  then  the  source  values  entered  via 
other  commands  will  be  used.  The  normalization  procedure  is 
described  in  Section  4.1.1. 

Examples : 

*N-NORM  2+23 

2 3 

which  normalizes  the  neutron  source  weights  to  2x10 
*G-N0RM  1 

which  normalizes  the  total  prompt  gamma  ray  source  spect- 
rum to  one. 

4. 1.2. 9 *y-YIELD  value 

y = same  as  in  Section  4. 1.2. 8. 

value  = weapon  source  yield  in  KT.  The  default  value 
of  the  yield  is  one  and  there  is  no  explicit  unit  definition  for 
this  command. 

Examples : 

*N-YIELD  100 

which  specifies  100  KT  as  the  neutron  yield  of  the  weapon. 
*G-YIELD  .02 

which  specifies  .02  KT  as  the  yield  for  prompt  gamma  rays. 

4.2  GEOMETRY  SPECIFICATION 

4.2.1  General  Geometry  Input  Description 

The  geometry  configuration  for  TDATR  is  illustrated  in 
Fig.  14  where  the  component  coordinates  are: 
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Hg  (HS)  = source  altitude 

(HT)  = target  altitude 

Rg  (RS)  = slant  range 

R^  (RH)  = horizontal  range 

9 (AN)  = slant  angle 

Three  consistent  component  specifications  define  a complete  | 

geometry  configuration  with  respect  to  the  ground  level  (as 
long  as  one  of  them  is  either  HS  or  HT) , and  the  other  two 
components  can  be  calculated  from  the  three.  TDATR,  therefore, 
requires  the  specification  of  three  geometry  components  one  of 
which  may  have  several  different  values  up  to  50.  The  specified 
output  results  are  then  displayed  for  the  various  geometry  com- 
binations. The  format  of  the  input  commands  is 

*xx,  units  , values 

where  xx  is  replaced  by  one  of  HS , HT , RS , RH , AN.  The  unit 
options  include  most  of  the  reasonable  units  that  are  appropriate 
(see  Section  4. 2.2.1).  Values  may  be  specified  in  a list 
separated  by  blanks  or  in  the  format:  n^^  (n)  n2 , signifying 
values  ranging  from  n^  to  n2  in  steps  of  n.  In  this  case  n2 
must  be  arithmetically  greater  than  n^.  Another  construct  for 
entering  values  is  of  the  form  r*v  where  r is  a repetition 
factor  and  v is  the  value  to  be  repeated  r times.  Any  of 
the  constructs  or  mixture  of  constructs  may  be  used  for  the 
definition  of  the  values.  The  description  of  these  two  con- 
venient constructs  is  most  appropriate  here;  however,  they  can 
be  used  with  any  other  command  that  requires  a list  of  values 
for  input . 

Two  of  the  possible  geometry  configurations  result  in 
ambiguities.  When  RH,  RS  and  HT  are  specified,  there  is  no  ; 

inherent  information  whether  HS  should  be  placed  above  or  ! 

below  HT.  In  order  to  resolve  the  ambiguity  the  characters 
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"+"  and  should  be  used  with  the  *HT  comniand  to  indicate 
that  HS  should  be  placed  above  or  below  HT  respectively  (e.g., 

*HS-  1000).  The  other  ambiguity  occurs  when  RH,  RS , and  HS  are 
specified  and  the  placement  of  the  target  height  is  in  question. 

In  this  case,  the  + and  - characters  are  used  with  the  *HS 
command  to  place  the  target  above  or  below  the  source  respectively. 
If  no  character  is  specified  then  the  default  is  +. 

There  is  also  an  option  in  TDATR  to  move  the  ground  to  a 
desired  altitude.  It  is  effected  by  the  *GR0UND  command  and  all 
specifications  involving  HS  and  HT  are  interpreted  relative  to 
the  ground.  When  sequential  runs  are  computed  with  TDATR,  the 
ground  is  not  automatically  reset  to  zero  unless  respecified 
by  another  *GR01JND  command  or  a *ST0P  command  is  encountered. 

4.2.2  Specific  Geometry  Commands 

4. 2. 2.1  *xx,  units,  value (s) 

XX  = one  of  RH,  RS , HT , HS  and  AN  denoting  horizontal 
range,  slant  range,  target  height,  source  height,  and  slant 
angle  respectively  (Fig. 14). 

units  = one  of  M (meters), KM,  MILE,  YD,  KFT , and  FT  for 
the  distance  specification  and  one  of  DEG  (degrees) , RAD 
(radians),  and  COS  (cosine)  for  the  angle  specification.  Be- 
cause of  the  symmetry  of  the  geometry  configuration  about  the 
source,  the  range  of  admissible  angles  is  from  -90°  to  +90°. 

Default  unit  for  the  distance  specification  is  M (meters)  and 
for  angle  it  is  DEG  (degrees) . 

value (s)  = one  or  more  values  of  the  corresponding  geo- 
metry component.  The  maximum  number  of  values  for  any  coordinate 
is  aroitvarily  fixed  at  50. 
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As  mentioned  above,  three  of  the  geometr>  commands  define 
a complete  TDATR  geometry.  Each  geometry  component  may  have 
several  values  in  principle;  however,  the  intended  use  is  for 
two  of  the  components  to  have  single  values  and  the  third  com- 
ponent to  have  one  or  more  running  values.  Another  typical  use 
is  to  have  multiple  values  for  two  components  and  one  value  for 
the  third,  in  which  case  the  two  components  with  the  multiple 
values  will  be  successively  paired  and  used  with  the  single  value 
of  the  third  component. 

There  is  a special  unit  definition  option  (GM)  for  xx  = RS 

which  is  allowed  if  the  following  geometry  configuration  is 

specified:  HS , HT,  RS.  In  this  case  the  values  associated  with 

2 

the  XX  = RS  command  are  interpreted  as  g/cm  of  the  slant  range. 
The  GM  unit  option  is  restricted  to  this  configuration  only. 

Examples : 

1.  *HT  +,  KM,  1 

*RS,  KFT,  .5  (.5)  10 
*RH  100 

which  specifies  the  target  height  at  1 kilometer, 
several  values  of  the  slant  range  starting  at  .5 
kilofeet  and  ending  at  10  kilofeet  in  steps  of  .5 
kilofeet,  and  the  horizontal  range  at  100  meters. 

Note  that  the  different  coordinates  may  be  specified 
in  different  lonits.  Also,  the  source  altitude  will 
be  calculated  to  be  above  the  target  altitude. 

2.  *HS  2000 

*AN,  DEG,  10  (10)  80 
*RS,  KM,  1 

which  specifies  the  source  altitude  at  2000  meters, 
the  slant  angle  ranging  from  10  degrees  to  80  de- 
grees in  steps  of  10  degrees,  and  a slant  range  of 
1 kilometer.  The  DEG  unit  specification  for  the 
slant  angle  is  superfluous  since  it  is  the  default 
unit . 

3.  *HS,  KM,  2 
*HT,  KM,  2 

*RS,  GM,  50  120  380  AlO  550 
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which  specifies  co-altitude  of  both  source  and  target 
at  2 kilometers  above  the  ground  and  the  slant  range 
of  5 different  values  in  units  of  g/cm2. 

4. 2. 2. 2 ^GROUND,  units  value 

units  = same  as  distance  specifications  in  4. 2.2.1. 

value  = distance  of  the  ground  relative  to  sea  level 
for  the  TDATR  problem  geometry. 

If  this  command  is  not  specified  then  sea  level  will  be 
used  for  the  ground  level.  If  this  command  is  specified, 
corresponding  HS  or  HT  specifications  are  interpreted  relative 
to  the  ground  level. 

Examples: 

1.  *GROUND,  KMf  . 5 

which  specifies  the  ground  to  be  at  a half  kilo- 
meter. 

2 . -'GROUND  500 

which  has  the  same  effect  as  Example  (1)  since  the 
default  unit  definition  is  meters. 

4. 3 OUTPUT  RELATED  SPECIFICATIONS 

4.3.1  General  Output  Command  Description 

These  commands  allow  the  user  to  specify  the  output  con- 
figuration of  results  from  TDATR.  The  *RESP  command  specifies 
the  type  of  response  that  is  required  for  output.  It  can  specify 
any  or  all  of  the  following: 


1. 

Total  Fluence 

2. 

Tissue  Dose 

3. 

Concrete 

Dose 

4. 

Air  Dose 

5. 

Silicon 

Dose 

6. 

Tantalum 

Dose 
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The  response  functions  used  to  generate  the  doses  are  listed 
in  Appendix  A. 

The  output  times  can  be  specified  either  as  local  time 
by  the  use  of  the  *LTIME  command  or  as  absolute  time  by  the 
use  of  the  *ATIME  command.  Only  one  should  be  specified  but 
both  appear  on  the  output. 

The  *TITLE  command  is  used  to  provide  the  title  for  the 
individual  TDATR  problem  and  it  is  useful  when  several  problems 
are  run  in  a given  TDATR  session. 

There  are  three  control  commands  that  are  also  described; 
*EXC,  *STOP , -'FIN.  After  the  source,  geometry  and  output 
commands  are  specifiea  for  a TDATR  problem,  the  *EXC  command 
is  used  to  execute  the  problem  and  display  the  results.  For 
subsequent  problems  only  those  elements  of  the  problem  description 
are  needed  to  be  specified  that  differ  from  the  previous  problem. 
If  radical  changes  of  problem  descriptions  occur  then  the  *STOP 
command  should  be  used  to  reinitialize  the  parameters  of  TDATR 
and  the  problem  must  be  described  in  full  detail.  The  *FIN 
command  is  used  to  terminate  a TDATR  session. 

4.3.2  Specific  Output  Commands 

4.  3.2.1  *RESP/z/(ij^  i,  . . .) 

z = either  G or  NG  representing  prompt  gamma  rays  or 
neutron  generated  secondary  gamma  rays  respectively. 

i^,i2,.-.  = response  option  indices  separated  by 
blanks.  These  are  single  digit  numbers  (1  through  6)  cor- 
responding to  the  fluence  and  dose  options  described  in  Section 
4.3.1.  Any  or  all  of  the  six  options  may  be  present  in  any 
order. 
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Examples ; 

1.  *resP/G/(1  2) 

which  specifies  the  output  of  total  fluence  and 
tissue  dose  for  prompt  gamma  rays. 

2.  *RESP/NG/(3  4) 

which  specifies  the  output  of  concrete  dose  and  air 
dose  for  neutron  generated  secondary  gamma  rays. 

4. 3. 2. 2 *LTIME/z/,  units,  values 
2 = same  as  in  4. 3. 2.1. 

units  = one  of  NANO.  SHAKE,  MICRO,  MILLI , or  SEC  repre- 
senting the  possible  time  units.  The  default  is  NANO  when  the 
unit  definition  is  not  specified  in  which  case  the  delimiting 
commas  must  be  absent. 

values  = local  time  values  from  low  to  high  order  re- 
presenting the  local  time  boundary  values  used  for  the  time 
dependent  output. 

Examples : 

1.  *LTIME/G/  .05  .1  1 10  25  35 

40  60  80  100  300  500  800 

which  specifies  the  output  local  time  boundary  steps 
for  prompt  gamma  rays  in  nanoseconds  which  is  the 
default  unit  specification. 

2.  *LTIME/NG/. SEC, 5-10  1-9  5-9  1-8  5-8 

1-7  6-7  1-6  3-6  6-6  S-6 

which  specifies  the  output  local  time  boundary  steps 
for  secondary  gamma  rays  in  seconds. 

4 . 3 . 2 . 3 ”ATIME/z/,  units,  values 
z = same  as  in  4.3.2. 1. 

units  = same  as  in  4. 3. 2. 2 including  the  default  para- 

mete  r . 
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values  = absolute  time  values  from  low  to  high  order 
representing  the  absolute  time  boundary  values  used  for  the 
time  dependent  output. 

Only  one  of  *LTIME  or  *ATIME  should  be  specified  for  a 
given  TDATR  problem  per  "particle". 

Examples : 

1.  *ATIME/G/ .SEC, 1-7  1-6  1-5  1-^  3-4 

6-4  103  102  .1  .3  .5  1 

which  specifies  the  output  absolute  time  boundary 
steps  in  units  of  seconds  for  prompt  gamma  rays. 

2.  *ATIME/NG/ .MILLI,  1-4  1-3  1-2  .1  .3 

.6  1 10  100  300  500  1000 

which  snecifies  the  same  values  for  output  absolute 
time  as  in  Example  1,  only  the  unit  specification  is 
different  and  it  applies  to  secondary  gamma  rays. 

4. 3. 2. 4 *TITLE  n 


n = up  to  74  characters  used  as  a problem  description 
title  to  identify  the  output.  The  title  will  remain  in  effect 
until  replaced  by  a new  "TITLE  command  or  cleared  by  the  effects 
of  the  *ST0P  command. 


Example : 

1.  '-TITLE  SAMPLE  TDATR  PROBLEM  1 

which  would  result  in  the  string  "SAMPLE  TDATR 
PROBLEM  1"  being  displayed  with  the  output  of 
TDATR  results. 

4.  3. 2. 5 --'-EXC 

The  effect  of  this  command  is  to  execute  the  TDATR 
problem  that  is  specified.  The  relative  order  of  the  previously 
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described  commands  is  immaterial  but  when  this  command  is  en- 
countered the  specified  problem  is  executed. 

4. 3.2.6  *STOP 

This  command  clears  all  flags  and  internal  buffer  areas 
set  up  by  previous  specifications.  TDATR  is  initialized  and  a 
complete  problem  specification  must  follow  this  command. 

4. 3 2.7  ^FIN 

This  command  terminates  the  program  with  a FORTRAN  STOP. 


55 


5 . REFERENCES 


1.  L.  Huszar,  L.  F.  Nesseler  and  W.  A.  Woolson,  "User's 
Guide  to  Version  2 of  ATR  (Air  Transport  of  Radiation) 
SAI-73-534-LJ,  DNA3144Z,  April  1973. 


2.  L.  Huszar  and  W.  A.  Woolson,  "Version  3 of  ATR  (Air 
Transport  of  Radiation),"  SAI-75-585-LJ , DNA3819F, 
July  1975. 


3.  L.  Huszar,  W.  A.  Woolson  and  E.  A.  Straker,  "Version  4 
of  ATR  (Air  Transport  of  Radiation),"  SA1-76-561-LJ , 
January  1976. 


4.  E.  A.  Straker,  W.  H.  Scott,  Jr.  and  N.  R.  Byrn,  "The 
MORSE  Code  with  Combinatorial  Geometry,"  DNA2860T, 
May  1972. 


5.  S.  A.  Dupree,  H.  A.  Sandmeier,  G.  E.  Hansen,  W.  W.  Engle,  Jr. 
and  F.  R.  Mynatt , "Time  Dependent  Neutron  and  Photon  Trans- 
port Calculations  Using  the  Method  of  Discrete  Ordinates," 
SAI-70-123,  August  1970. 


6.  W.  H.  Scott,  Jr.  and  B.  L.  Colborn,  "FLAIR  - A Scaling 
and  Folding  Code  for  the  Generation  of  Photon  Transport 
Results  in  Air,"  SAI- 74- 516-LJ , March  1974. 


APPENDIX  A 


TABLES  OF  PERTINENT  PARAMETERS 


Table  A-1.  Source  and  detector  energy  boundaries  for 
prompt  gamma  rays , and  detector  energy 
boundaries  for  secondary  gamma  rays  (MeV) . 


1. 

0.02 

- 0.05 

10. 

1.33  - 

1.66 

2. 

0.05 

- 0.1 

11. 

1.66  - 

2.0 

3. 

0.1 

- 0.2 

12. 

2.0  - 

2.5 

4. 

0.2 

- 0.3 

13. 

2.5  - 

3.0 

5. 

0.3 

- 0.4 

14. 

3.0  - 

4.0 

6. 

0.4 

- 0.6 

15. 

4.0  - 

5.0 

7. 

0.6 

- 0.8 

16. 

5.0  - 

6.5 

8. 

0.8 

- 1.0 

17. 

6.5  - 

8.0 

9. 

1.0 

- 1.33 

18. 

8.0  - 

10.0 




Table  A-2.  Source  spectra  for  neutrons. 


Group 

Energy  Boundaries 
(MeV) 

Weapon  Fission 
Source  Values 
(Fraction 
in  Group) 

Thermonuclear 
Source  Values 
(Fraction 
in  Group) 

1 

1.07(-7)- 

2.9(-5) 

0 

0 

2 

2.9(5)  - 

1.01(-4) 

0 

2.00(-3) 

3 

1.01(-4)- 

5.83(-4) 

0 

2.40(-2) 

4 

5.83(-4)- 

3.35(-3) 

0 

1.22(-1) 

5 

3.35(-3)- 

0.111 

2.227(-l) 

3.65(-l) 

6 

0.111 

0.55 

1.693(-1) 

1.02(-1) 

7 

0.55 

1.11 

2.159(-1) 

8.50(-2) 

8 

1.11 

1.83 

1.468(-1) 

6.20(-2) 

9 

1.83 

2.35 

1.060(-1) 

2.80(-2) 

10 

2.35 

2.46 

5.743(-3) 

5.00(-3) 

11 

2.46 

3.01 

2.871(-2) 

1.90(-2) 

12 

3.01 

4.07 

5.481(-2) 

2.60(-2) 

13 

4.07 

4.97 

1.177(-2) 

1.70(-2) 

14 

4.97 

6.36 

1. 832(-2) 

1.80(-2) 

15 

6.36 

8.19 

1.274(-2) 

1.47(-2) 

16 

8.19 

10.0 

7.342(-3) 

1.41(-2) 

17 

10.0 

12.2 

0.0 

2.56(-2) 

18 

12.2 

15.0 

0.0 

7.06(-2) 

[ 

i 
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Table  A-3.  Dose  response  functions  for  gamma  rays  ( rad/ ( y/cm  ) ) . 


F 


Z-N 

Z-N 

Z-N 

Z-N 

Z— N 

Z-N 

Z-N 

z-s 

E 

z-s 

z^ 

z-s 

o 

o 

o 

o 

O 

o 

o 

O 

z-s 

z-s 

z-s 

z— s 

z-s 

z^ 

p 

ON 

ON 

ON 

ON 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

I-H 

C^x 

a> 

On 

eyv 

CJN 

<7N 

1— 1 0) 

1 

1 

cd  cn 

S-z 

S-z 

V-z 

N«Z 

N*Z 

'wZ 

N-Z 

N-z 

s-z 

NwZ 

NwZ 

NwZ 

NwZ 

s-z 

SwZ 

SwZ 

S-Z 

SwZ 

u o 

CN 

r*^ 

rH 

rH 

o 

CO 

CO 

UO 

CO 

CO 

>d- 

rH 

CM 

CTv 

d- 

C Q 

00 

m 

UO 

NO 

UO 

00 

•d- 

rs 

o 

o 

CO 

CO 

CO 

CO 

. 

• 

» 

. 

. 

. 

. 

• 

. 

H 

oc 

NO 

CO 

rH 

ON 

NO 

UO 

UO 

UO 

NO 

CO 

rH 

rH 

CM 

CM 

CO 

UO 

Z-N 

z^ 

Z— N 

z^ 

Z-N 

Z— N 

z— N 

Z-N 

z-s 

Z-N 

C 

o 

rH 

rH 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

z-s 

Z-N 

z-s 

Z— N 

z-s 

o cd 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

ON 

C5N 

ON 

ON 

C^ 

o B 

1 

1 

S*z 

v-z 

s-z 

NwZ 

N-Z 

NwZ 

NwZ 

s-z 

S-z 

NwZ 

NwZ 

SwZ 

s-z 

S-z 

SwZ 

SwZ 

rH  a» 

ro 

LO 

UO 

UO 

00 

o 

UO 

O 

CN 

O 

UO 

o 

00 

CO 

00 

o 

•ft 

rH 

CM 

rH 

VO 

CO 

CM 

rH 

o 

rH 

rH 

00 

CM 

st 

00 

CM 

oo 

C/D 

. 

. 

• 

. 

. 

. 

• 

nO" 

ON 

rH 

rH 

CM 

CO 

UO 

NO 

C- 

CO 

ON 

rH 

I-H 

I-H 

CM 

CM 

z-v 

/~s 

Z-N 

Z-N 

Z-N 

Z-N 

Z-N 

Z-N 

Z-N 

Z-N 

Z-N 

z-s 

rH 

rH 

o 

o 

O 

o 

O 

o 

o 

o 

O 

o 

Z— N 

z-s 

Z—N 

Z-N 

Z-N 

CO 

f-H 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

I-H 

ON 

CTn 

C3N 

ON 

CJN 

1 

1 

•H  C 

N-Z 

S-z 

v—z 

VwZ 

N-Z 

N«Z 

NwZ 

NwZ 

NwZ 

NwZ 

S-Z 

SwZ 

NwZ 

NwZ 

s-z 

NwZ 

SwZ 

S-z 

< <u 

o 

CSJ 

CM 

rH 

CO 

00 

NO 

o 

UO 

UO 

CO 

o 

CM 

CM 

o 

O 

d- 

be: 

COv 

rH 

vO 

CO 

CN 

rH 

O 

rH 

rH 

CO 

UO 

I-H 

CO 

vO 

<3N 

CM 

vj- 

<N 

to 

rH 

rH 

CN 

CO 

UO 

NO 

CO 

I-H 

rH 

rH 

rH 

CM 

Z-N, 

z^ 

Z-N 

Z-N 

z^ 

Z-N 

Z— N 

Z— N 

Z— N 

z-s 

z-s 

0) 

o 

O 

rH 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

z— s 

z-s 

z-s 

Z-S 

z-s 

iJ 

r-1 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

rH 

I-H 

rH 

rH 

On 

CJN 

ON 

(On 

CTN 

o td 

1 

1 

1 

1 

1 

I 

(-1  B 

N-Z 

S-z 

V«-z 

\-z 

s-z 

SwZ 

NwZ 

NwZ 

NwZ 

s-z 

SwZ 

s-z 

s-z 

SwZ 

SwZ 

NwZ 

s-z 

O i-l 

o 

O 

O 

O 

00 

CN 

o 

o 

UO 

UO 

UO 

O 

o 

00 

VO 

O 

00 

UO 

c u 

<7N 

CN 

o 

CM 

v£> 

<t 

CO 

rH 

o 

rH 

I-H 

vJ* 

00 

I-H 

Nt 

00 

rH 

vD 

o ^ 

• 

. 

. 

. 

. 

. 

. 

o 

iri 

rH 

00 

rH 

rH 

CN 

CO 

UO 

NO 

CO 

CJ^ 

I-H 

rH 

I-H 

CM 

CM 

c 

z~s 

z^ 

Z-N 

Z— N 

Z-N 

z— s 

Z-N 

Z-N 

Z— N 

Z— N 

z-s 

o 

rH 

rH 

rH 

O 

o 

O 

O 

o 

o 

o 

o 

o 

z-s 

z-s 

z-s 

z-s 

Z-N 

z-s 

CO  0) 

t—i 

rH 

rH 

rH 

i-H 

I-H 

rH 

rH 

rH 

rH 

I-H 

rH 

(T' 

ON 

ON 

ON 

CJN 

CJN 

M p 

1 

1 

1 

» 

(U  CO 

N-Z 

NwZ 

N.^ 

N-Z 

N_Z 

N^ 

NwZ 

NwZ 

NwZ 

s-z 

SwZ 

N-Z 

NwZ 

SwZ 

SwZ 

s-z 

SwZ 

SwZ 

•a  M 

o 

O 

CM 

vO 

O 

UO 

o 

UO 

UO 

00 

Cn. 

ON 

VO 

CM 

C -r^ 

cn 

ON 

v£> 

CM 

uO 

UO 

CO 

CO 

o 

CM 

UO 

O 

<f 

d)  H 

< 

. 

. 

• 

. 

. 

. 

00 

m 

VO 

rH 

rH 

CN 

CO 

nJ- 

UO 

NO 

OO 

rH 

rH 

I-H 

rH 

CM 

CM 

>^ 

M bOz-N 

m 

CO 

NO 

(U  Vj  > 

o 

rH 

CM 

CO 

NO 

00 

o 

CO 

VO 

o 

UO 

o 

o 

o 

UO 

O 

O 

Q,  0)  (U 

. 

. 

. 

• 

. 

. 

> 

. 

. 

a c S 

o 

o 

O 

o 

o 

o 

o 

rH 

rH 

rH 

CM 

CM 

CO 

UO 

VO 

00 

O 

DU'-' 

rH 

a 

3 

rH 

CM 

CO 

UO 

NO 

00 

CTN 

o 

I-H 

CM 

CO 

UO 

NO 

Cn. 

00 

O 

rH 

rH 

I-H 

rH 

rH 

rH 

rH 

rH 

rH 

Vj 

CJ 

i 

1 


59 


Wl 


Table  A-4.  Detector  local  time  boundary  values 
for  prompt  gamma  rays. 


Group 

Time  (sec) 

Group 

Time  (sec) 

1 

0.000 

14 

1. 000(-7) 

2 

1.000(-9) 

15 

1.468(-7) 

3 

1.468(-9) 

16 

2. 154(-7) 

4 

2.154(-9) 

17 

3 '2(-7) 

5 

3.162(-9) 

18 

4.642(-7) 

6 

4.642(-9) 

19 

6.813(-7) 

7 

6.813(-9) 

20 

1. 000(-6) 

8 

1.000(-8) 

21 

1.468(-6) 

9 

1.468(-8) 

22 

2. 154(-6) 

10 

2.154(-8) 

23 

3.162(-6) 

11 

3.162(-8) 

24 

4.642(-6) 

12 

4.642(-8) 

25 

6. 813(-6) 

13 

6.813(-8) 

26 

1.000(-5) 
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Table  A-5.  Detector  absolute  time  boundary  values 
for  secondary  gamma  rays. 


Group 

Time  (sec) 

Group 

Time  (sec) 

1 

0.0 

22 

1.78(-4) 

2 

2.15(-7) 

23 

2.05(-4) 

3 

4.64(-7) 

24 

2.37(-4) 

h 

1.00(-6) 

25 

2.74(-4) 

5 

1. 78(-6) 

26 

3. 16(-4) 

6 

3.16(-6) 

27 

4.22(-4) 

7 

5.62(-6) 

28 

5.62(-4) 

8 

1.00(-5) 

29 

7.50(-4) 

9 

1. 33(-5) 

30 

1.00(-3) 

10 

1.78(-5) 

31 

1.58(-3) 

11 

2.37(-5) 

32 

2.51(-3) 

12 

3.16(-5) 

33 

3.98(-3) 

13 

4.22(-5) 

34 

1.00(-2) 

14 

5.62(-5) 

35 

1.78(-2) 

15 

6.49(-5) 

36 

3.16(-2) 

16 

7.50(-5) 

37 

5.62(-2) 

17 

8.65(-5) 

38 

0.10 

18 

1.00(-4) 

39 

0.316 

19 

1.15(-4) 

40 

1.00 

20 

1. 33(-4) 

41 

2.00 

21 

1.54(-4) 
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APPENDIX  B 


DENSITY  SCALING  OF  THE  TIME  DEPENDENT  BOLTZMANN 
TRANSPORT  EQUATION  (BTE)  FOR  NEUTRAL  PARTICLES 


In  this  appendix  we  prove  the  general  scaling  law  for  the  flux 
that  exists  between  systems  having  material  densities  which  are 
distinct  but  which  are  constant  in  space  and  time.  Specifically, 
let  the  flux  in  system  1,  (r^  , t , E , n)  , be  known,  and  let  it  be 

desired  to  calculate  the  flux  in  system  2,  <t>2  (^2  ’ ’"Z  ’ ^ 
the  following  relationships  are  true 


P2  = 

r2  - k r^ 
t2  = k-^t^ 


J 3 — 


j 3 — 


S2(r2.t2.E.fi)d  r2dt2  = S^(rj^,t^,E,Q)d  r^dt^, 


'^2^^2’^2’^’^^  = k $j^(rj^,t^,E,n)  . 


then  we  conclude 

c 4>j^(r 

To  show  this  we  write  down  the  BTE  for  system  1 
^ 9<J'^(r^,t^,E,fi) 


3t, 


(1) 

(2) 

(3) 

(4) 

(5) 


S2^(r2^,t^,E,G) 


Q‘V^<t’^(r^,t^,E,fi) 


- Pj^-aj.(rj^,E)  •<I'^(r^,t^,E,n) 


+ p^-/dfi'/dE'aj.(rj^,E')  -GCE'  ,n'  ,E,n)  . 


'fl(ri,ti,E'  ,n) 

where  the  notation  is  conventional. 


(6) 
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mCEDINa  PUX  Burac-NOT  FILMBD 


We  now  make  the  substitution  indicated  by  Eqs.  (1-4)  into  Eq.  (6). 
Note  the  following: 


9 _ 1 9 

9t,  k 9t, 


"l  = k "2 


Also  Eq.  4 can  be  written  as 


Si  (Ti  , ti , E , = ~^S2  (r2  , t2  , E , 

k 


Thus  we  have 


V^k^9t2‘*’l^^l’’^l’^’^  ^ 


— 2(r2>t2tE,n) 


- (i)0-V24>^(r^,t^,E,n) 


i P2  • (r2  .E)  • , t^  ,E  ,fl) 


i P2-/dfi'/dE'aj.(r^,E') 


• G(E'  ,51!'  ,E,J^)<J>[,r^,tj^,E'  (10) 


Now  if  we  make  the  substitution 


^ ~ k^$j^  (r^  , t ,E  ,5^) 


into  Eq.  (10),  we  get  the  BTE  for  system  2 and  our  conclusion 
follows . 


APPENDIX  C 


The  time  dependent  radiation  environment  is  determined 
by  a convolution  of  the  time  dependence  of  the  source  and  the 
time  dependence  of  the  radiation  transport  results  obtained 
for  a delta  function  source.  That  is 


where 


oCt) 


S(t)  {>  (t-i) 


di 


S(t)  is  the  time  dependence  of  the  source 
0 (t-i:)  is  the  time  dependence  of  the  transport  of  radiation 

and 

J)  (t)  is  the  desired  radiation  environment. 


The  following  figures  represent  some  plausible  examples  of  the 
functions  involved  in  the  convolution.  The  functions  S(t)  and 
(|)(t)  represent  the  source  function  and  the  transport  function 
respectively  in  the  first  two  figures.  The  third  and  fourth 
figures  represent  the  pictorial  view  of  the  components  of  the 
convolution.  The  third  figure  shows  the  transport  function 
(i)(-T)  as  rotated  about  the  axis  and  the  fourth  figure  shows 
the  same  function  shifted  by  the  constant  t^  and  superimposed 
on  the  source  function. 
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given  t = t 

° f(T)  = S(t)  0 (t^-T) 


Although  both  S(t)  and  (()(t)  are  treated  as  histograms 
in  the  code,  they  can  be  thought  of  as  continuous  functions 
by  properly  interpolating  on  the  histogram  values. 


Typically,  S(t)  is  a smooth  continuous  function  such  as 
a Gaussian,  and  this  analysis  is  based  on  that  assumption.  The 
nature  of  ^(t)  is  such  that  its  time  resolution,  relative  to 
S(t)  may  be  very  coarse.  Therefore,  if  we  consider  the  range 
of  integration  as  being  over  S(t),  for  theoretical  considerations 
we  can  take  <t)(t)  to  be  slowly  varying  over  the  range  of  S(t)  and 
thus  a method  of  integration  considering  functions  which  are  no 
more  complicated  than  S(t)  is  adequate. 

Experience  has  shown  that  for  most  smoothly  varying,  well- 
behaved  functions  the  method  of  gauss  quadratures  performs  very 
well.  The  12-point  Gaussian  quadrature  can  be  written  as: 


6 

I 

i=l 


w. 


a + c(l+x 


i>]  " 'I 


a + c(l-x 


i>]| 


(1) 


66 


where 


= interval  weights,  0 < x^<  1 

6 

w.  = amplitude  weights,  I w.  = 1 

i=l  ^ 


The  essence  of  the  problem  is  the  representation  of  the 
various  functions  involved.  Normally,  all  three  of  the  functions; 
S(t),  <))  (t)  , ^>(t)  have  different  time  resolutions  and  therefore 
the  interpolation  method  may  be  important.  In  order  to  preserve 
the  maximum  overall  structure  a three-point  Lagrangian  inter- 
polation is  generally  adequate.  In  some  cases,  however,  a two- 
point  (i.e.  linear)  interpolation  may  be  adequate  but  for  further 
accuracy  a higher  order  spline  interpolation  may  be  required.  A 
parabolic  (three  point)  interpolation  is  an  adequate  compromise  for 
most  of  the  problems  being  considered  and  is  thus  utilized  in  TDATR. 

If  we  let  t^  represent  the  discrete  time  values  and  let 
y^  represent  the  corresponding  function  values  of  S(t^)  or  <})(t.) 
then  the  form  of  the  three-point  Lagrangian  interpolation  is 
as  follows:  given  points:  (t^,y^),  and  a 

point  t such  that  t^  < t < t^  and  |t-t2l  £ min  (It-t^l, 

|t-t^|}  where  |••.1  denotes  the  distance  norm,  then  y 
corresponding  to  t is  given  by: 

(t-t2)(t-t2)  (t-t^)(t-t2)  (t-t^)(t-t2) 

^ ~ ^1  (t^-t2)  (tj^-t^)  ^2  (t2-tj^)  (t2-t2)  ^3  (t^-tj^)  (!t3-t2)  ‘ 

This  scheme  forms  the  parabola  defined  by  the  three  given 
points  and  evaluates  the  point  y corresponding  to  the  given  t 
from  the  parabola.  Since  the  solution  of  the  convolution  integral 
is  an  important  part  of  TDATR,  it  is  useful  to  consider  the  amount 
of  computer  time  consumed  by  this  procedure. 


The  following  components  of  the  total  time  must  be 
taken  into  account; 

T^;  time  of  interpolation 

T2;  time  of  functional  evaluation 

T^:  time  of  integration 

The  interpolation  time  can  be  minimized  by  keeping  track  of  the 
interval  and  if  it  is  the  same  as  before  then  certain  inter- 
mediate variables  can  be  saved  and  need  not  be  computed  every 
time.  There  are  some  other  overall  efficiency  considerations 
which  could  be  exercised  but  will  not  be  discussed  further. 

If  the  floating  point  multiply  of  a computer  is  used  as 
a tinit  of  operation,  an  add/subtract  as  a half  unit  and  a 
divide  as  four  units  then  an  estimate  of  the  computation  time 
can  be  obtained  in  units  of  multiplication  times.  The  follow- 
ing estimates  for  the  three  times  are  obtained  by  simply 
counting  the  number  of  operations  in  the  above  equation  and  in 
Eq,  (1): 


T^  = -^  (add/subtract)  + 9 (multiply)  + 3*4  (divide) 

= 28  units 

which,  in  principle,  has  to  be  evaluated  for  every  point. 

T2  = ^ (add/subtract)  + 4 (multiply)  + 4T^ 

= 118  units 

this  gives  the  number  of  units  necessary  to  evaluate  the 
quantity  in  Eq.  (1)  which  is  in  brackets. 
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= 7 (multiply)  + 6T2 

= 715  units. 

There  are  other  overhead  operations  which  will  increase  the 
time  somewhat  but,  on  the  other  hand,  if  there  are  functional 
or  parametric  representations  of  S(t)  and  (p  (c)  then  the 
interpolation  time  can  be  greatly  reduced. 

For  the  Univac  1108  this  analysis  gives  about  1.9  msec 
per  one  point  of  '5>(t).  If  the  interpolation  is  performed  in 
log  space  then  an  extra  20  units  must  be  added  to  for  the 

exponential  function.  This  gives  a total  of  1195  units  or 
about  3.1  msec  on  the  Univac  1108  computer.  This  would  entail 
that  both  S(t)  and  p (t)  exist  in  log  space  as  well  and 
that  logs  need  not  be  taken  at  every  step  of  the  interpolation. 
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The  following  seven  sample  problems  constitute  a repre- 
sentative set  of  the  capabilities  of  TDATR.  The  first  three 
examples  are  concerned  with  the  secondary  gamma  rays  due  to 
different  neutron  sources.  The  rest  of  the  problems  show 


examples  of  prompt  gamma  rays  due  to  sources  of  various  con- 
figuration. 


Sample  Problem  1 


This  problem  presents  the  results  due  to  a neutron 
fission  source  for  secondary  gamma  rays.  Two  responses: 
tissue  dose  and  silicon  dose  are  presented  at  a geometry  con- 
figuration of  source  and  target  heights  at  1 km  and  a separ- 
ation of  500  m.  The  output  absolute  times  are  entered  in 
units  of  seconds. 
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TDATR  SAMPLE  PROBLEM  1 - INPUT 


*TITLF  f TS.STilN  SnilHCF  StCC'MUPY  r,AMMA  RAYS  Af  IKM 
♦ N-SdL'RCF  1 0 0 
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B.65-S  i.is-'i  i.s^-a  «?-'!  2,7n-u  a,pp-a  7,'^-ii  t.ss-s  a-3  i.h-2  S.^-a 

.1  .-Mb  1 l.RRR 
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TDATR  SAMPLE  PROBLEM  1 - OUTPUT 


rc'41‘<  Prt(.)BLt;M  MOviHfW 


MSSIiIN  SDOPCf  SKrOMOAWY  gamma  rays  at  1km 


*»****ll*l,t***********»»************ 


(gf-MTRUN  SMHRGE  iRM  At  I / A 1 I ( l^=  I.OdOE-fOn  mMI!K('NS/kT 

Yltin=  i.onor*”"  ''i  tiuai  iiiiiput=  i .onoL  + oo  ui.utRONS 

REUiMis  - hissii'N=  uodiiiKno  Thh P'-ni‘'iiiri.E ARi  o.ono(;-oi  wimiv=  o,noot-ni 


* * * * 

*.  \ 

* * * * Tk 

*t*»**»tt 

* ♦ * ♦ * 

♦ 

TDA  1 R Pki|B|  F M 

F I SS  I()M 

Sfl'iwrl  sFrilMBAKY  GAmmA 

WAYS 

* « • * 

r.Riutuo 

* * * 

1 1 VFI 

* * 

•*«*«* 

11  , tllloK  t1  , 

(1  .nnor.M/tP*  *?, 

»»«»-* 

n . npoKF  I , 

*<*•** 

n . 0 (1 0 M 1 1. 1 S 

A 

*HllR  1 /.  . 

pamgf 

WM  = 

(l.'v)fiKK', 

'>'i  .SB  FC.m/Cm  **?, 

1 .MlOKF  1 , 

n . 1 11 M I L E S 

SI  ANI 

KANI'.F 

HSr 

0,‘.1iiriKM, 

SS.sn^r.M/cM*  *?, 

1 .mHOKFT, 

n . ? n M 1 L F S 

T A(.(.F  t 

AFT. 

M - 

1 . 0 '1 0 » M , 

1 IM. 

S.PBI‘<F  I, 

<1.631  MILFS 

SOIImC  F 

A)  T. 

1 ,ih1(ikF>. 

1 in.  /01'',m/Im«*?, 

A.RB 1 kF  T , 

n .63 1 M I LE  s 

*SLAN I 

AM-I  f 

A K'2 

0 , fiiiniH,(,iu  F5  (C'is=  i.onoon) 

• (Atn'lAlFII  FRI.M  DlUfP  milRO  ? NA  rts 


4 

4 4 4 * « 4 * 

4 4 4 4 4 

4 4 ft  4 4 ft 

* ft  4 4 ft  4 

4 ft  ft  4 4 4 

4 4 ft 

.';F  (.nSliARY  GAmmA  1 I SSI'f 

nnsF  (RAO/sFr) 

ahmu  nTF 

1 IT.  Al. 

T I “F 

1 I MF. 

iincF’EI.  K'F  0 

Scat  If  PFii 

tiiIaL 

ClIMMLATIVt 

} 

<. ToOf -06 

3.0  53f  "II6 

0 . nonc-o 1 

1 .6ftnt-1 7 

1 .6F10L-1  7 

6.  l«6t-23 

4 

O.MSoF  “(16 

/.  763E -06 

o.nnnt-nt 

1 .S loF -t  7 

1 .Sl'lE-l  / 

t .«20t-33 

s 

I .BSnF 'oS 

1 .66  <F  -cS 

0 . OonF  - (H 

6.01 iT -|B 

6.0  1 1E-1FI 

3.aoSF-32 

b 

F.3'»ST  -OS 

A.  1 36F -oS 

0 . oooE -0 1 

«. F3nE -iR 

o.saoE-iR 

3.SRRt-33 

1 

S.  <6iiE-oS 

S.  1 U -r,S 

n • 0 0 0 E - 0 1 

'F.6Q|E-3n 

'>.oo,E-3n 

?.f>22E  -22 

fl 

7.S/SF  -ftS 

1 , - OS 

n . 0 0 n F - 0 1 

O.RPlF-pO 

O.R61F-30 

2.633E-22 

q 

1 . 00  /r -ou 

^I'S 

o.onnF-0( 

6. 1/OF -30 

6.370F-3O 

2.6S0F-33 

10 

1 . WSE-O'J 

1 . <36E -OU 

n,oooF.-n  1 

6.0S0F. -30 

6.0SOF -30 

2.67'»l-33 

1 1 

l./ZoF-na 

1 . 7S  SF  -no 

n . onnE-n 1 

6.0601 -30 

6.  nORF -30 

3.  /n2E-33 

13 

P.  WOF  -0  1 

3.  <S  U -0« 

0 .nnnF_-o  \ 

6.06AF -3O 

6.0(,  01.-30 

3.70  7f -33 

I A 

A.uRnF  -n<4 

< . 06  iF  -liO 

n . oonF -0 1 

S. 7o  U -30 

S.  /<4  IF  -20 

2. "321-33 

1<A 

S.F'o  af  -no 

o.oooE-01 

S.237E-30 

S.33Tt-30 

1.00  lE  -33 

I . ( 6SF  -n  S 

1 . 1 6 <F  -n  A 

0 . 0 0 n F - n 1 

n .6  /6F -30 

0 .6 /6F -30 

3. SR|E -33 

IM 

P.TOnt-ok 

3 . Tl'MF  -n  A 

0 . nooF -0 1 

0 . n /6E  -30 

0.0  /HF_-3n 

0.3/HE-33 

1? 

1 . I noE-o3 

1 . t ilOF  -03 

0 . 0 0 n F - n 1 

i.n7HE-3n 

l.o7Ht-30 

F1.6PHE  -33 

1« 

S.  TpoF -03 

TnnF -02 

n . onoE -0 1 

1 . GSRE -30 

1 .6SRF:-30 

I .0001 -31 

I** 

7 .B(, (iF  - f)3 

A.«noF  -113 

n . oooE -n  1 

6.630E -31 

6.630E -31 

|.7Rit-21 

?ll 

P.OF^oF-Ol 

P.nBoF -01 

o.nooE-ot 

I . i7nE-31 

1 .3/OF-31 

2.0671 -31 

31 

6 .SBoE -0 1 

6 .SMoF -0  I 

n.oooF-Ol 

S.031t-35 

S.023E-33 

?.12|E-31 

?3 

1 .SooE  +00 

1 .oRiTF  Knn 

o.oont-ni 

t .63RE-32 

1 .63FIE-32 

?.2"«t-3i 

« 

• * # # 4 4 4 

4 ft  4 4 4 

4 ft  4 4 4 4 

444444 

444 
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TDATR  SAMPLE  PROBLEM  1 - OUTPUT  (Cont'd) 


StCOMOABY  GAmmA  SIlI''-0n  OOSfc  (HAD/StC) 


AnsnuinE  I fJCAi. 


T I^E 

t 

iiNcnLLintr) 

SCATIEHtn 

TOIaL 

CUMULATIVE 

i 

^ “06 

^».o^a)--o6 

O.OOOE-OI 

?.?0SE”17 

^.aoSE-l 7 

8.160E”^} 

a 

‘J.asoF-oh 

7.  7M^E-06 

n.OOOF-Ol 

1 .6«/l-l  7 

1 .6R7E-I  7 

P.56flE-22 

s 

1 .«5oF-oS 

1 .hH  if  -oS 

n.nonf -oi 

7.6«71-lH 

7.6fl7E-jn 

T.  168E-PP 

b 

-05 

i.  I?8I:-0S 

o,0ftOE-oi 

1 .0((hE-l» 

1 ,Oa6F-t0 

1. 561E-2P 

f 

S. ShOf -05 

S.tOH-OS 

o.ooor-ol 

1 .PP^F-l'l 

1 .22iF-i6 

3.  189E-22 

H 

7.S?5F-05 

7,(4()«F-o5 

o.nnoF ”01 

S.OO  u -P0 

5.<7a^E-?n 

3.on?E-22 

Q 

1 .007F  -0'< 

O.Oiifit  -PS 

o.onot-oi 

7.01  iF-r’O 

7,a|  lF”PO 

3.74PJL-PP 

10 

1 .^asr-oo 

1 . ^?f*F-00 

0 . 0 n 0 1 - 0 1 

6. ^n7i -PO 

6.  ^o7E-?p 

5.oa7E-22 

1 1 

1 . 7 / 0 >•  - 0 't 

t . 7S  -oo 

0 . 0 0 0 1-  - 0 1 

6,  lO'Jf  -PO 

6.  to'7E-?n 

3.077F-PP 

?.  57ot-  -no 

P.  1,'t  tt-iiu 

0 • 0 0 0 1 ”01 

6.  \pot-pn 

6.^P'7F”P0 

3.SpoE-PP 

t < 

} . oMoH -00 

06  U -(i.'i 

0 . 0 0 n 1 ”01 

S.60  -PO 

S.'»oiE-po 

3.7m  pt  ”22 

n 

S,«noF-nO 

S . 0 0 i I-  - 0 u 

0 . on  pi -0  1 

S,0S71 -po 

S.0S7E-PP 

3.  7'7it-2P 

1^ 

1 . 1 65K -o  5 

1 . t -(1  \ 

0 , oopf-  - 0 1 

a .flHsi -po 

a ,RRSt -po 

0.  197E-P2 

1 

70oF-ni 

p.  70«t -ni 

o.oopr-nt 

0 .p6  5t-po 

O.P6lF-pa 

5.PprtL-pP 

1? 

1 . 1 OOF -0^ 

t . 1 onf -0^ 

0 . 0 0 (1 1 - 0 1 

l.PPPl”PO 

i.?P^F”PO 

9.7  \Q1 -p? 

l8 

<./0c|F-0? 

^ . 7i)oE-{iP 

0 . il  0 .1 1 ■ 0 t 

1 . 7 -PO 

1 . 7 ^OF-po 

1 .6  ^at -PI 

7.0nor-o<» 

7 . n n 0 ^ - 0 p 

0 . pool ”0  1 

h.os/iF-pi 

s.osaF -pi 

1 .9aol -PI 

^0 

^ . 0»^nf  -1, 1 

^ 1 

0 , pp  pf  - 0 1 

1 .0  -PI 

1 .0  iOF-pi 

p.p^lt -PI 

6 . 0 fi  - 0 1 

n.S'^iiV  -n  1 

0 , 0 ,1 1 - 0 1 

S.p/fl7 -p  1 

s.p7;'f  -p  1, 

P.P07t -PI 

?.<> 

1 , S 0 n r.  ♦ 0 0 

1 .00  ((  *(10 

o.nopi -pt 

1 . 7l pF-pp 

1 . 7 1 pF -pp 

p.aSRt-pi 

* r X?.  c'i  T I UN 

* * * « * 

F !f  l> 

*****  ft 

****** 

* * * 

Sample  Problem  2 


This  problem  is  the  same  as  the  first  one  except  that 
it  presents  the  results  from  a thermonuclear  neutron  source. 
Since  only  the  source  and  title  changed  from  Problem  1,  those 
are  the  only  input  specifications  necessary. 
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TDATR  SAMPLE  PROBLEM  2 - INPUT 


*^ITL^  THlWMIlMin  f- AW  S Uiwrf  S^r.O^^OAWY  r,AMMA  WAVS  AT  1KM 
■*N-snuKCF  0 t n 
*1  vr 


ru*-i ^ 


TDATR  SAMPLE  PROBLEM  2 - OUTPUT 

! 


TPATW  PWdHl.fM  NUMlitP  2 THfWHONUCI.EAW  SOURCfc  StCONO&BY  IJAHMA  RAYS  A)  IKM 

4**4**«*A-k*****4*4A«***fti(*AA*A****A 

NfiiiRnN  siiiiRC'-  niiR'iAi  ] /AT inN=  i,ooin,»no  nkni rons/k t 

YIKin:  1.0i10(-4(10  Kt  riHAL  lUIPlMi  l.oonk+on  ^ltl)TW(^NS 

wt  iMtrS  - FISSI'W-  n.OOflk-RI  iHtRK'Ml'Cl  k ARr  |.0(int+0n  IRMkVr  O.POOE-O) 


*********A«*A******4«******«******* 

IPATW  PRflUifM  HIImBFR  2 1 MFWMIlKlHri  k AR  SnURkk  SkCnMOARY  r.A.,hA  RAYS  AI  |km 

«4lt*«««**ft»*A*A*A*4«*A4k**lkA**A*AA* 

pRiiiiAn  I f VH  rt.iioor.M/r.M»*?,  n.nnnKFT,  o.ono'UkkS 

•Howl/,  RAhlPr  WMr  U.SnA^M,  SS.SmGR/kM*  <?,  t.kUl'KFt,  11  . i I 1 M I |.  F S 

SI  AMI  KANl.F  Rk':  i1,S1>PKm,  sS  . SB  SG’a/Cm*  • , l.pnOKFI,  B.MIMILES 

TARSI  I alt.  tn  - t.OlinKM,  n h.  7RI  r.M/oi*  *?,  a.phiimi,  ii.I'Pimilfs 

sni'pri  A(  r.  hs=  i.oor.xM,  n b.  7‘T  i hm/Cw  ^,^«UkT,  o.spihTlIS 

• f;.  amt  amf.i  f am-  a . nciPi'i'i.Pi  k s (nisr  t.onOPii) 

‘FAI.CUI  A IFI;  FROM  mi'FR  kill  IWII I M A T k S 

ft*«*ft*«*««***«il*ft******ft**t«**-***** 


Sir.M'glMWY'  flSSl'^  ()f’5»'  (MAO/SCC) 


AHSfU  lll^  \ Iirai 


T 1 Mk 

r } Mf- 

iRiriHl  infii 

SCM  TFREII 

T i'll  A L 

CIIMIII.  A n Vk. 

J 

7rtnF-oS 

£ • 0 -OK 

n . oonk  -0 1 

1 . A'l/k  -Ih 

1 . 7->/r-lb 

b.«2/k-P2 

a 

R . 'jSflk  -(>M 

f • ^^3^-  ’~()h 

0 • onnk -0 1 

1 .20  At  - |b 

I .20  Al  -Ib 

I .BoRl -PI 

s 

1 , 0 S '1  k * n S 

1 . O -nS 

n , noi'k  "n  1 

0 . 0 1 '!  F - I A 

j . n R '1  k - 1 7 

,?  . 0 iSk  -pt 

kl 

^.;)RSk  -liS 

A.  IP'^k  -CiS 

n • 0 n n k - 0 1 

<1.  WBE-l" 

'l.  A7Bk  - jH 

p.  t 1 bl  -PI 

1 

S , i k»  (1  k - n S 

s.  1 0 AI  -ns 

n . nnnk  -n  t 

A.I  .A7E-1R 

A.I  A7k-lR 

2.  I2AE-P1 

fl 

k.s/sf  'OS 

A . iOBk  -Ii'i 

n . onnk  - n t 

B . 'JR  5k  -?n 

B..IR  Ak  - 20 

2.1PSk-pi 

R 

1 . no  'k  -T/J 

O.RnBi  -oS 

i>. nook  -0 1 

b . I JRk  - po 

b . 1 JRL-PO 

2. IPbk -PI 

1 0 

1 . k'lSk  -O'l 

1 . ApBk  -f,U 

0 . nook -0 1 

S.BB|k-pO 

S.BHtc-PO 

2.t2Rk-Pl 

I 1 

1 . / /nk -nR 

t . AS  AI  - no 

n . oook-o 1 

S.7i7t-P0 

S.7i  Ak-pn 

2.1 Alk-Pl 

1 2 

p.  A/nk-o'i 

p.  AS  Al  -on 

0 .non'  -0 1 

s,  7 ( sk  -po 

s . A ( Sk -pn 

p.  1 ABl  -PI 

1 ^ 

A.  iBoF-oU 

1.  lih  Ak  -O'J 

0 . nook -0 1 

S.'isnk  -po 

S , U S 0 k - p 0 

P.  1 'J'lk  -PI 

1 '1 

S.Ohflk -n i 

s . « ) Ak  -n'l 

0 . oonk  -0 1 

S.oPAE-po 

s.opAk  -Pn 

2.1''0k-pi 

IS 

I . 1 o'lk  -nS 

) . 1 '■  il  -0  A 

0 . n 0 0 1 - 0 1 

J.SARk -pn 

•J.SA'Ik  -20 

2.1 RBk -PI 

1 k* 

p.  iRok-n  A 

P.  ABBk  - fl  < 

o.oook-oi 

A.OSBk -PO 

A.'ISBk  -PO 

2.2R  Ik  -21 

1 . 1 noF-  p 

1 . 1 Bl)k  -OP 

0 . onnk  -0 1 

p .'oRBk  -pn 

p.RROI  -Pn 

2.  Al  Ak  -PI 

1'’ 

A . Aook -OP 

A.  7. ink  -op 

0 . ionk-o I 

1 .MR  AI  -po 

1 .SR  Ak  -PO 

A.  AiRt  -21 

1'’ 

/ . Book  -oP 

A .Bunk -op 

0 . oonE*0 1 

b.  Ap  Al  -PI 

b.  Ap  Ak  -Pt 

A.SR7k  -PI 

20 

p . oBok  -111 

p.nBiik  -ni 

n.OOoE-pi 

t . A0Bf-P1 

1 . AoBk-Pt 

A.O /Of -21 

21 

B , S B 0 k - 0 1 

B . S B n f - m 

0 . nook -0 1 

J . 7 a j k - p A 

II . Au  jk  -P3 

A.'MPk-Pt 

22 

1 .SOrtk  K10 

1 . 'IR  Ik  1 0(1 

O.OOoE-01 

I .SpBF-pp 

1 .SpOk -pp 

a . nb JE-Pt 
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TDATR  SAMPLE  PROBLEM  2 - OUTPUT  (Cont’d) 


SECnNOAwy  C.AmmA  SltlCfUi  DOSF  (WAH/SPC) 


AMsn|.n^^  inrAL 


I I«E 

iiNcHLUntn 

SCiFIEHFO 

tqTal 

CIIMMLArlVE 

5.70OF-nh 

O.OOOF-Ot 

1 .fla'iF  -16 

1 .TlaPF-lO 

6.95ot-aa 

a 

‘>.'4bfiF-06 

7.7e^F-Oh 

O.rtOCF-OI 

1 .aspF - 1 p 

l.aSPF  -16 

1 .660E-at 

s 

1 .""jni  -0*’ 

1 "OS 

O.OnoE-ol 

«. ia7F-| 7 

a.5a7F.-l  7 

a.iipfc-ai 

s 

^.^‘>SF-oS 

i.  1 -o'* 

p.dOoF ”01 

«.63at-tP 

a.63aE-i H 

a.i'76E-ai 

7 

S. ?60E-OS 

S.  1'41F-0S 

n.oooF-ot 

i.a«6F-ip 

i.ap6F  -19 

a.aoit-ai 

8 

7.S7SE-05 

7 . Ud^E-oS 

(1 . npiiF-o  I 

6. EppF -ao 

6. 7fl0F -ao 

?.apsF-ai 

9 

1 .no7F-oa 

g.g„HF-oS 

O.OOPF'Ol 

p.aiaF -ao 

6.?iaF -ao 

a.ao6t-ai 

to 

t . i«SF-0'l 

1 . sa^E-o" 

0 . OOpF  "P  I 

S,Ph7F-ao 

5.967F-ap 

a.ap9F-ai 

1 1 

I . 7 /nF  -n'l 

1 . 7S  <F  -o'l 

p . oont -0 1 

6.PP6F.-ao 

f>  tOohL-pn 

a.ai it -ai 

\? 

s7nF-0<i 

a.  u.-fu 

0 , aooF-*’  1 

h.pppt-ap 

6.Pn6F-aO 

a,ai6E-ai 

I i 

; . uHoF  -nu 

A.'ihif  -0*i 

0.00  OF -PI 

s,  7 $1 F -ao 

s. 7}lF-ap 

a.aaoF -ai 

1 

^ • ^hot-  ~rni 

s.rt'l  U -,t'4 

P . IIOPF  -0  1 

s.aESF -ao 

s.a«st  -ao 

a.aoaF -ai 

i'- 

1 . 1 hSE-n  i 

1 . 1 MF  -Hi 

0 . OPOF  -0  I 

P . 7HPF -aO 

« . 7«oF -ao 

a.asiF -ai 

ix 

?.  7'>0F  -.1^ 

a.  7PHF  -(>5 

P.POOF  -01 

■1.17  iF  -ao 

C 

EV 

1 

a. isaF -ai 

1 r 

) . t OflF  -(V 

1 , 1 nn(  -0^ 

0 , 00(U  -0  1 

i.i6SF-ap 

i.l65F-aO 

a.oasF -ai 

i« 

70(lF  -iV 

/oot -n^ 

P .PPOF -p 1 

1 .6«'IF-?0 

1 .6Kat-ap 

i.'iesf  - ai 

/ , "0^* 

/.‘'mF-oa 

P . PPPF -P 1 

p.ppaF-ai 

p.poaF-ai 

1.76oF-ai 

<»o 

/> , 0 n h - (*»  1 

a.nanF -ni 

p . POOF -0 1 

! . iHPF -ai 

1 . iP6F -ai 

p.osoF-ai 

^1 

^ . S'*nF  -0  1 

P . PPPF -0 1 

s.p  i'>F.-a  1 

s.p  1F.F  -ai 

'i.popF  -ai 

1 . S (1  n t ♦ 1)  0 

1 . .'i'F'FF  Ki(\ 

P .pool  -P  1 

1 .6a  iE.-aa 

1 .6a  iF.-aa 

u.as6E-ai 

* * 

****** 

****** 

****** 

****** 

* * * 

**  r n luN  roMPi  ► Tf  o 
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Sample  Problem  3 


This  problem  is  the  same  as  the  first  one  except 
it  presents  results  from  a 12.2-15  MeV  neutron  source, 
only  the  source  and  title  specifications  are  necessary 


that 

thus 
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TDATR  SAMPLE  PROBLEM  3 


INPUT 


TITLf"  I^-MTV  S^r(I^JDA^<Y  fjAMrA 

0 I)  1 

F xr 


WAVS  AT  Ihm 
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■ ■ <4 


TDATR  SAMPLE  PROBLEM  3 - OUTPUT 


TDATW  PWOHLfH  NDmhKR 


I'J-MEV  SFCCIMOARY  gamma  mays  at  iKM 


*********^*k*****k**itt************* 

nkiiiron  sniiijfF  mhwm  al  ' t ln'M=  t.oooF+on  nHjirons/kt 

YIFir)=  l.OPOFtOn  Kl  iriTAL  m.lIPUTs  i,opof>oo  nH'TWONS 

«tI(;HlS  - F1SSII1N=  O.nnpF-01  THt  KM()NUCLFAR=  O.OnOF-oi  liiMtVs  1.00uE*00 


* A « A 

AAA 

A A A 

A A A A 

»*»<,**** 

***** 

A 

A A A A A A 

A 

TDATR  PWiTRi  Fm 

NUmHI-  p 

4 

l 'J-Mt  V 

SFCdMDAkY 

gamma  pays  at 

1 KM 

• * » * 

* * « 

A A A 

A A A A 

******** 

***** 

****** 

A 

GRI'DM) 

1 i-vru 

0 

. OOOKM  , 

0 . 0 (1  n (>  M ^ C M * * p , 

IT . 0 n 111'  F 1 

p.nAfiMii  ts 

*Hi)P  t ^ . 

wA^.^.^ 

0 

,Sl»OKM  , 

SS .SM  TGm/Cm*  *p , 

1 . M ii  0 " F T 

<1 . 5 U M I IJ  S 

SI  AN  I 

a'am;f 

s = .) 

^ -■)  0 K M , 

SS.SkSHR/C^^aP, 

1 .MUOKF  t 

1) . i 1 1 M I L F S 

I A Pl.t,  r 

Al.  1 . 

1 

, 0 (>OK  M , 

M M . / y t G M / r , 

T .^P 1 KF  1 

n.hPlMIl FS 

Si  il'Wl  t 

4|  T . 

HSr  1 

,n*»OKSi, 

1 ) M . / a 1 c.M/r  y * ♦p. 

T,P8  1*F  t 

IT  .G<T|  mile  S 

* S 1 A I'l  1 

A NCI  F 

A\  = 

II . iTpoijF  r.Pf  F s ((;i)S=  i.noono) 

A * « « 


AfAicuiArf-n  oiMf-H  cofiwo  t NA  u s 

s'-f'fiNi'iwY  t.amma  iissuf-  nosi-  (w^n/stc) 


AMSOLiiiE 

1 IlFM- 

TTMF 

r I F'F 

ilflGHI.I.  IDF  D 

Sc  A I 'F.FFFD 

kitaL  r.dMniAilVf 

5 

S • TOPf -os 

P . P 5Pf -ps 

n . 0 Pill  -p  1 

1 .S«5F'|5 

1 .SA  Af  - )S 

S.  A'FSt  -PI 

'J 

s . IIS|)(  -II  s 

7 . 7PpF -OS 

n , PiipF  - II 1 

«. /PoF-lS 

A.  7SoF--)S 

1 . AlsF -ao 

S 

1 . PSn F - ii S 

j.hPiF-fiS 

n . p p II 1 - 0 1 

7>.  5'I'iE-iS 

,1 . A 'i  a r - ] 6 

1 .SSpt  -pci 

6 

5.?RSF -ns 

t pPF  ~0S 

n . PiipF  - p 1 

af-i  7 

P.Pl  AF  -1  7 

t .Sl^t -^0 

7 

S. TSpF -nS 

S . T P S 1 - 1 1 S 

n . p 0 p F - 0 1 

p.  iS'lF  -1« 

IS  IF  - iS 

1 ,Sl 7F -^n 

n 

/ .S  /',F  -OS 

7 . II  pPF.  -liS 

n . p 0 II F - n t 

1 . /7PE-1S 

1 . 77VK-,P 

1 .SI 7E-P0 

'} 

t . 0 0 / r - (I  u 

P.SliPF  -pS 

n , p p p F - n 1 

p . AspF -^0 

?.  ASnt  -rip 

t .SlAL-^0 

1 0 

1 . VJSF  -O'T 

t . sppf  -pp 

n . p n p F - p 1 

1 ,'7p  AF;-,'P 

1 . '>Pi^-PD 

I .SlflF -20 

1 1 

1 . 7 / 0 F - 0 0 

1 . 7 S S F - (1  U 

n . npot -PI 

1 .rtP7E-^’0 

1 . AP7F  ~i>r, 

1 .St AF-^n 

IP 

P . W 0 F - 0 '■* 

P . 3S5F  -P'l 

0 . p 0 p F - n ( 

1 .t'SlF  -^0 

T .AS  It -^0 

1 .St  Ar-20 

1 < 

S.yPnf  -o'l 

T. JS  U -PU 

n . nopF -0 1 

1 . 7 1 At -^0 

1.7] AF -^0 

1 . S 1 A 1.  - 2 P 

I'J 

S.Phiil  -Oj 

S.Pil  U - C/4 

0 . POOF -P 1 

1 .S 1 SF -^0 

1 .StSF-,>0 

t .S 1 RF -^0 

IS 

1 . 1 SSI -0  < 

t . t »'  At  - 0 5 

0 . nopF  - ii  1 

1 . Ap^t  -,'0 

1 . AiiMt  -^0 

1 .S20L-2n 

1 *' 

P.  /''pF  -(li 

P.  7PPF-0  5 

p . iippF  -0 1 

1 . 1iT1F-?0 

1 . IPlF-^0 

1 .S22F -20 

1 ' 

1 . ,0Pr-P? 

, . , pp|;-p^ 

p . 0 p p F - n 1 

«.  ISMF  -.'1 

A.  ISAF -?1 

t .S5'tF  -20 

IH 

T , 7 P n F - 0 ^ 

A . 7 n i|F  - r,  p 

p . pppi -p 1 

p. ASiF -pt 

a . 5slF -^1 

1 .ssoF-2n 

/.PpoF”l)»* 

7 . PPiiF  -lip 

n . POpF -p 1 

1 . /'HF  1 

1 . 7.1AI  1 

1 .SSAE-?o 

?0 

p . n p p F - p 1 

p.pppt  -fll 

n . piipF  - 0 1 

A.S5iE-p^ 

A.hAlF.-^? 

1 .SSSE-2P 

^1 

S.S«of  - III 

s . S “ oF  - II ) 

n . p p 0 F-  - 0 1 

1 . AASF-<>5 

1 . AASF-^5 

1 .Ss7t -20 

P7> 

1 . SqoF  * n p 

1 .'I'F'FF  *,,0 

0 . OPPE -0 1 

'1  .P'j  af  -pi 

a . aa  At  -Pi 

1 .S/lF -20 

TDATR  SAMPLE  PROBLEM  3 - OUTPUT  (Cont'd) 


Sf-COUCARY  Hll,  fCONi  DUSE  (RAIl/StD 


AUSHLiiTE  I.OCM. 


TIMf.- 

t IMF 

iiscm.i.  tof  D 

Sr/F  > Ih  WK) 

TflTAt 

ClIFHILA  I I VE 

3 

3 . 7ooE 'Oh 

2.0i2F -Oh 

0 . OOOF -0 1 

1 • /21F-|S 

1 . /21E-1S 

h.S39L-21 

a 

R . uSnE 'OR 

7.  7B^F  -pt, 

0 . 00  of "0 1 

<>.  52hF  -ih 

R.32«’F-16 

1 . 3721 -20 

s 

l.flSilt-o5 

t .hF-  5t-05 

0 . OOof "0 1 

2.SuoE-1 h 

2.SooE-th 

1 .636t-20 

S 

3.2R‘>E-o‘3 

3.  1 2hE"0S 

O-OOOF-Ol 

3.07SF-1 r 

3.07Sf-i7 

1 .h95E-20 

7 

S , 3M17  -OS 

S.IRU-OS 

n,oooF-oi 

2. 3S'|F  -1<3 

2. 3SnF -18 

1.6997-20 

ft 

7 .S /Sf -OS 

7.'|nHF, -OS 

0 , OOOF  -0  1 

1 .Ry  3F -in 

1 .R53F-1O 

1 .6997 -20 

1.O07F-0J 

R.OOXt -OS 

o.onoK-01 

2.S2'’t  -20 

2.S2Rt-20 

1 . 6997 -20 

in 

I . -on 

I . 32''/  -O'l 

0 . 0 0 0 E - 0 1 

2.0  5'>F  -2n 

2.0  3hE-2n 

1 .699L-20 

1 1 

1 . 7 /OF  -O'l 

1 ,7SU-0u 

n.n(ioF-0| 

2 , OoHF -20 

2.0087-20 

1 .6997 -20 

12 

2 . Wol  -O'l 

2 . 3S  3(  -0/4 

0 . 0 ri  0 F.  - 0 1 

t .‘77oF -20 

1 .R70F -20 

1 .6991  -20 

1 3 

3.<JI*oF  -o'l 

f . uo  3F -on 

0 . 0 0 0 F " 0 1 

1 .rt| OF -20 

, .H,8F -20 

1 . 7oo7 -20 

1 ■* 

s.«00(  -n.’i 

S . «'/  FL-O/i 

0 . OOOF  -0 1 

1 .h02F -20 

1 .h02F.-2n 

1 . 7O0t -20 

1 '•> 

1 . 1 OSF -0  \ 

1 . 1 h 31,-0  5 

0 . 00  OF -0 1 

1 . 3 7hF -20 

1 .5/hF-20 

1 .7017-20 

1 

2.  /‘’oF  -0  S 

2.  /'F’iF  -0  3 

0 . OOoF  -01 

1 . 1 h 3F  -20 

1 . 1 h 3F  -20 

1 . /O'lF  -20 

1 > 

1 . 1 0 0 F - 0 2 

1 . 1 'oF  -02 

0 . OOoF  -0 1 

8 . hr.OF  -2  1 

8 . so  8F -2 1 

1 . 7 1 oF -20 

7ooF-n2 

3 . /ooF -02 

0 . OO'iF  -<t  I 

'1.  SO  /F.  -2  1 

'1 . SO  /F  -2  1 

1 . 7.3  JF  -2  0 

r’ 

7 , « 0 O F - 0 2 

7,0007  - 0 2 

0. oool -0 1 

1 .rt.l2'  -21 

1 .8.12F  -21 

1 . / '421  -20 

20 

2 . 0 n 0 F - ,1 1 

2 . 0 H 0 F - 0 1 

0 . 0 0 0 F - 0 1 

3.H22F -22 

3.822F -22 

1 . /SoF -20 

21 

. SOoF  -0  1 

'1  . S.3.)l  -0  1 

A , oonF -01 

1 . '10  'll  “2  3 

1 .OO'IF  -2  3 

1 . /’>!  1 -20 

22 

1 . F. 00  F ■»  0 0 

t . /I'lRF  * on 

0 . OOoF  -0 1 

■■1  .hhOf  -2  3 

'1  . S'. ‘IF  -2  3 

1 . /'.6l  -20 

« 

* ^ 

**««•«« 

xh  ruT  Tn».  rtivpj 

• 4 « « * 

»■  I rn 

* * 4 4 * 4 

4 ft  4 
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Sample  Problem  A 


This  example  presents  results  due  to  an  8-10  MeV  prompt 
gamma  ray  delta  function  source.  The  requested  responses  are 
the  concrete  and  air  dose  at  1 km  co-altitude  with  a source- 
target  separation  of  800  m.  The  output  local  time  mesh  is 
given  in  seconds,  and  it  will  be  used  for  all  subsequent 
problems . 
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TDATR  SAMPLE  PROBLEM  4 - INPUT 


♦ TlTLf-  H-in  MFV  sniJRtr  PWfjwpT  gamma  rays 

*HS,KM, 1 
*HT,KH, 1 
ARS  800 
AG-SilUWCF  ( 1 ) 
aG-SVAL  1 MO  1 

*LTIMP /G/,SFr.,  0 t-R  a.hU-9  |-h  ^.is-H  a.(S4i-fl  1-7  2,lS-7  <|  6u-/ 

t-6  «?.1S-6  ii,ba-h  1-b  ‘ ' 

♦ RESP/G/C^  <4) 

AEXC 


TDATR  SAMPLE  PROBLEM  4 - OUTPUT 

TD»1K  ►'POHLFM  NUMBfW  I n-io  MFV  SOIJRCF  PwnwPT  GAHMA  PAYS 

GAHmA  SnUPTF  MIPM  A|  I z A I I (I'jz  l.nOOE  + OO  GAhmAS/kI 

'^IFinr  1 .000f«ii0  KI  niTAL  OUtPIlIz  l.noOF*00  gammas 

FIMt(SFr)  / ENFi?(,Y  1 FNFPGY  e tMFPGY  < FNtPGY  U (NFPGY  5 FNFPGY  6 
O.ooooF-01  n, (loot -01  o.o(.iu-oi  ij.ooni-(ii  (i.O(.ot-oi  o.iionf-ot  o.onoL-Ol 

IIMFtSFn  / 1NFW1.Y  / f^EPlJY  B F'IFKI.Y  tlJtWItY  | 0 FGEPGV  t(  fNFPGY  \^ 

O.(ino0(.-ot  o.iinn(-o)  n.iiimf-oi  n.oiiot-oi  n.onni-nl  o.oont-ni  o.ooof-ot 

**tt*4*»*»t*««***A4*«****A4A**AA*A* 

TMtl'^Ff)  / FNFPI.Y  I a FNtPl.Y  1 'I  F>aW(,Y  FNtW(;Y  Ih  FNtPl.Y  J/  tNtPl.Y  (M 

o.noooF-oi  0.0 II (It- 0 1 n.iKiof-oi  n.nnot-ni  fi.ooof-oi  o.ooof“01  i.ooolaOO 


roAiw  (>w,ili|  t s.  ijiiMUt  u 1 H-to  ut  V sniii-rt  pkoimi  G'mmA  pays 


(.PIIUM) 

1 f VI-  1 

(1  . oco  - “ . 

0 . onoi.'vr*-’  • 

0 . 0 0 (1 » t 1 , 

n.oOpMii 1 S 

*“llP  1 / . 

>J  f.  f- 

0 . POCk "» 

AM  . o <p(,,4/r  .1  • 

* ^ » 

P . opSi^F  1 . 

(1 , aO  /m  I ( t s 

SLaM 

WS: 

.•.„or.i.M, 

H«,‘tAP(.M/rM* 

O.qPFMIt  f s 

I A Pl.t  t 

M T . 

H ] - 

1 .00(1''.  M, 

1 1 (>  . / 0 1 I.M  /C  M • 

n.MP|“U  F S 

S'  ll'F  < F 

H T. 

1 , (")(ln  » 

1 1 M . /o  t Gm/C^* 

^.PMUF  I , 

0 . M P 1 M 1 ( F s 

• S L A M 

ANbI  li 

A*)r 

O , I'lOrF  I.P(  f S ( ( dSr  1 . 

on  0 00) 

• r.M  CM.  A iki'  iiiHtP  r.iiiipn  I NA  If  s 


TDATR  SAMPLE  PROBLEM  4 - OUTPUT  (Cont'd) 


pri.impt  T.ahma  COMC'f^EtE  OdSl-  («AO/sfC) 


AUsni  ute  10C4L 


1 I 

1 1 MF 

HUf.liLL  InFh 

SCAI  HPFD 

iriTAi- 

CIIMlil  A 1 IVE 

1 

P . fetVf -nh 

s 

. OOOI  - 1 0 

8. 1 3Zf.-12 

n.opot-o  I 

5. 137F-12 

S.  1 57t-21 

p. 

. 8 7 pf-  - 06 

1 

.S/Zt-pP 

O.OOOF-01 

5.2297 -1 5 

3.2297-1 3 

8.5097-21 

i 

-08 

. Zf  -09 

O.OOOF-Ol 

2. 5067  -1  i 

2. 3P8F-1 3 

8.0037 -21 

a 

^,8787.-08 

7 

, 3201 -09 

O.OOflE-OI 

1 . 5u 17 -1  5 

1.5917-15 

8.8017-21 

s 

a.8''9l  -08 

1 

.S7S7-C« 

0 . OOoE-0 I 

8.8207-10 

8,629E-1  (1 

7.58  57  -21 

6 

P.70Pt -08 

3 

.59‘>I  -0« 

o.onoE-ni 

2.0887  -1  a 

2.0S8E-19 

0.2797-21 

7 

P. 7upf -08 

/ 

. 3207 -no 

0 . onoF-o 1 

1 . 1 5 1 7 - I 0 

1.1517-19 

0,8«0t-2l 

H 

2.O?8t-08 

1 

.S7St-o7 

o.ooot-ni 

a. 5227-15 

9. 5227-15 

9.3777.-21 

9 

5.n08K-08 

3 

-o7 

o.oooE-oi 

1 .8997-18 

1 .899E-18 

9.«Oo7-21 

1 0 

5.aniE-o8 

7 

. 5207 -0  7 

0.  OOOE-O 1 

0,8307.-18 

0,5307-18 

1 .0287-20 

1 t 

9 .?0'tF-n8 

1 

.S787-08 

o.ooo7;-oi 

9.0877-16 

9.HS7E-18 

1 .0827-20 

l^ 

8,o8'ir-o8 

3 

. 3987 -o8 

0.0007-01 

1 .0787  -1  8 

1 . 0 787- 1 8 

1 .1?"7-2o 

13 

‘),9«9K-p8 

7 

. 3201 -08 

0 . ooot -0 1 

1 .5287-1  7 

1.5287-17 

1 .1587-20 

ppuMPr  r.AHMA  A IP  niisi  (WAo/sf  C) 


AOSDI  i)lE 

I I f'7 

Ll'CAI. 

T I M(- 

IINC'U.L  tOF  n 

.5c  A 1 IF  PEP 

TOfAL 

Ct'OlPFATlVF 

1 

2 . 889j -08 

8.i)noi-  -10 

9 . 3 9 2f  - 1 2 

0 . r,  p p r - 0 1 

9. 592r-l 2 

9. 5927 -21 

2 

2 . n 7 0 1 - 0 8 

I .8Z77 -09 

0 .0007  -0 1 

2.821F-1 5 

2.0217 -1 3 

9.8877 -21 

5 

2,8  7,>^  -on 

3.  39  ZI  -(.9 

0 . npo7 -0  t 

2.0  387  -1  3 

2.0387 -1 3 

8. 1 Z 37 -21 

9 

2 . h 7 8V  -0^ 

7 , 1207  -09 

0.0007 -0| 

1 .205t-| 5 

1 .2087 -| 3 

5.8197-21 

8 

2 . 8 >'  J 3 - 0 8 

1 .6  ZSf  -08 

0 . 1 0 0 7 - 0 1 

8.0897-19 

8,08n7-ia 

8.5197-21 

h 

2. Zp9f -08 

S.  39', 1 -9O 

0 . O()0‘  -0  1 

2.89f,F-i  9 

2.8987 -| 9 

z. 1 9o7 -21 

Z 

2 . Z 427  -08 

Z. 3207 -08 

0 . OOOF - t 1 

I , P « 97  - 1 9 

1 .0891  -1  ;j 

Z. ZZar-21 

M 

2 . n 2 8 4 — p 

1 .8ZS7 -07 

0 . OOOf -01 

9.1981-18 

'1  . 1 9 87.-  1 8 

8.2517-21 

9 

3 . ni''‘7  - p8 

3 . 39S(  -It  Z 

0 . IH'Ol  - 0 1 

1 .8np7 - 1 8 

1 .8007 -18 

8.82  (7  -21 

1 0 

3 . 1 0 1 1 - 0 8 

'.  32'il  -4  7 

P . non! -0 1 

8 . 1 9 if  - 1 8 

8.19  37  -1  8 

8.9li  3'  -21 

1 1 

9 .2  197  -,18 

1 • Z )7  - 0 8 

0 . 0 0 0 7 - 0 1 

1 , 8fi87  — 1 6 

1 . 888F - 1 8 

9..19/|1  -2  I 

12 

8 . p89l  -n8 

3.  3981  -(18 

0 . 0(I07  -0  1 

'i.28nl  -1  z 

'4 , 2807.  - 1 Z 

9.2not -2I 

1 3 

9.9^91 -n6 

z . 5c’'’7  -08 

0 , oooi -0 1 

2.8  Z ZF- I 8 

2.8Z7E-18 

9.2157-21 

• • f » I r IP  r I n»!  rn'‘p|  i iM) 
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Sample  Problem  5 


This  example  is  similar  to  sample  problem  4 with  the 
only  difference  being  that  the  source  gamma  is  in  the  1-1.33 
MeV  source  bin,  and  only  the  source  and  title  specifications 
are  needed. 
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TDATR  SAMPLE  PROBLEM  5 - INPUT 


*TITLf-  1-1, :i3  M^v  SllJkCf-  PPHMPT  gamma  WAVS 

♦R-SVAL  fl*o  1 W*n 

*FXC 
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TDATR  SAMPLE  PROBLEM  5 - OUTPUT 


Trtirw  PHnHL^M  2 Mtv  souwct  prhmpt  camma  rays 

»*t****tttA»******»»t*****k******1l» 

r.AnnA  SOURCF  MIWmA)  t/A  r T(1N=  l.OOOE  + 00  GAMMAS/KT 

viti.ii=  i,oiiof+on  Kf  IDIAL  uutpijti  i.onnt'  + oo  gammas 
riMhtSFC)  / ENKWGY  1 CNiK(,Y  p KNKRGY  J energy  <J  energy  S energy  6 

n. onnoE-01  n.onoE-oi  <A,rion(_-oi  o.nonE-ni  n.nr.oE-ni  n.nooE-oi  o.oooe-pi 

*»***tt*»*ii********ttti,ti*i,t**t**** 

TIME(SEt)  / ENERGY  / ENERI.^Y  fl  ENERGY  q ENERGY  10  ENERGY  t|  ENERGY  IP 

o. ooo(iE-ni  n..ionE-oi  o.oooe-oi  i.nnoEfOO  o.oooe-oi  o.oooe-oI  o.oope-o| 

IINE(SEr,)  / ENEEJGY  li  energy  lu  EEIERGY  IS  energy  1G  ENERGY  17  ENERGY  IB 


O.OflOtiE-ol  0,il00E-01  O.OOOE-OI  O.OOOF-Ol  O.OOOF-OI  O.OOOF-01  O.OOOE-01 


* * « * 

ft  ft  ft 

ft  * 

ft  ft  ft  ft  ft 

ftftftftftftftftftft^ftftft* 

ft  ft  ft  ft  ft  ft  * 

rnAiR  pwnBi^FM 

NltM'EFR  P 

1 - 1 . 13  mf  V SOURCE 

PROMPl  gamma  rays 

ft  * « * 

r.PU"^'r'’ 

ft  ft  ft 

L^  VM. 

A ft 

ft  * ft  ft  ft 

0 . n 0 0 H , 

ftftftftftftftftftftftftft* 

0 , ooot.‘*/r^*  ft  n,oon*<^i’. 

ft  ft  ft  ft  ft  ft  * 

o.oon^lLtS 

/ . 

rAnIF 

e^M  = 

0 . rook  M » 

HB.R  lpi.M/CM**P,  P.BPSkFT, 

0.aR7MlLES 

SI.A'-I 

ANf^h 

O.B'IOnM, 

BB.‘E3Pr.N/rM**P,  P.aPSkfI, 

O.aq/MiLE  S 

TARI-l  T 

All. 

HT, 

1 .ono'<". 

116.  I.PPINFI, 

0.6P1MII tS 

SIEI'RE  f 

A|  T . 

"s  = 

1 ,(iaoKK» 

116. /RIGm/Cm**,*,  I.PMIkEI, 

0.6PtHlLf  s 

*Si  ant 

ANt.l  F 

A't  = 

0 .oonriE  GRE  F 5 (ri'S=  1 .000001 

• rALF'iiArFo  eriim  uiheh  cheeroinaifs 


* * * 
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TDATR  SAMPLE  PROBLEM  5 - OUTPUT  (Cont'd) 


PROl^PI  CONtPtir  r)iiSb  (PAO/StO 


APsnlDTP  LOCAL 


T 18p 

1 I 

iiMcOLL  inri) 

SCaTTFRFD 

tiiial 

cumlilat  1 VE 

1 

;>,8n9(.-n8 

S.OOOf -10 

5.  <0«F“iy 

O.OOOE-OI 

<0«E“1  a 

^.i08F“2i 

? 

^,8  /nL-08 

1 .57/L-oP 

o.OOOl “01 

1 .?iaL-iS 

1 .ai«b-is 

3.a«8F -23 

J 

^.87?P-o8 

J.  i0/f-(,9 

0 . 00oE“0 1 

?.1<>0E“1‘’ 

?.iaoE-is 

3.97St-2i 

u 

? . 8 78L  “()6 

7 . 5;?0t  “O'* 

0 . nooE“o 1 

1 . 666E- 1 S 

1 .666E-15 

u.868E-^3 

s 

-06 

) .S7Sb -08 

O.OOOK-Oi 

1 .4!78F-iS 

1 .278F-1S 

6.33aE“2i 

6 

7o^r-n6 

^.^6St-n8 

O.OOOt-Ol 

7.  1 7;>F  -16 

7.172F-I6 

8. 12aF-23 

1 

a.7a?F-06 

/.5^oE-0« 

0.fl00f“01 

6.0^^^ -1 6 

6.02iE-l 6 

1 .li5E-22 

8 

a.«?8b -0^ 

1 .S7SF-0  7 

0.O00F“Ol 

i.?i8t-ih 

5.?5«E-ib 

1 .So?E-aa 

9 

'S  . OOPb  “06 

5.^0SL-n7 

O.OOOF“01 

?.?«?E“l6 

2.?«2E-16 

2.076F-22 

10 

i.aoii  “O'* 

7.  5i!0F-07 

0.000fc“01 

2.07RE-16 

?.07<»F-i6 

3.1O0E-22 

1 1 

a.?oaL-o6 

1 .57SF-06 

0 . OOPL-O 1 

1 .Sa8F-l6 

1 .5il8l-l6 

a.970E-22 

IR 

6 . 06i|L  -ob> 

iObL -06 

O.OOoL-fll 

S . 0 i 1 E - 1 7 

S.0^lb-l7 

6.223b -22 

1 i 

R.pH9L-(,6 

7.  i£>0f  “06 

0 . 0 0 0 1 “ 0 1 

S . 11  1 E - 1 rt 

5. 1 1 lt-i8 

6.0971-22 

* 

• ***•* 

****** 

*•*»** 

****** 

****** 

* * * 

PNflMPT  GAh^'A  AIk'  OfiSt  (PAD/seC) 


AHRfUDlL  l.llCAL 


1 1'-t 

T !■'( 

"'ICMLI.  Kib  n 

SCAbTLWFn 

TOTaI. 

tllPIO-A  b 1 VE 

1 

2.689f-iU> 

, 0 OOF  - 1 (1 

3 . 30  8f  - 1 0 

0 . 0 0 n b - 0 1 

3 . 3nBb -1 0 

3.  3081  -2  3 

2 

2. 67 01 -06 

1 .S  7 />  -o» 

0 . 000b -0 1 

1 .21 2E-1 S 

1 .2l2b -iS 

3.008L-23 

5 

2. ''721-06 

3.  b'O/b  -o<' 

0 . OiipF  - 0 1 

2. 1 1 7b  -i>) 

2.  11  7b -)S 

3.9/0b-2S 

■'1 

2.878L  - ,t6 

7. 320b -pO 

0 . 0 0 0 b - 0 1 

1 . 6 h 1 F - 1 5 

1 .661*  -16 

f,.86SF-23 

S 

1 .6  /'>b  -08 

0 . 0 0 0 b - 0 1 

1 . 2 7 1 b - , S 

1 .27lb -1  S 

6.  326t -2  3 

6 

2.  /02b  -07^ 

3.  39'>b  -08 

0 . op  ob -0  1 

7 . 073t - 1 6 

7.0  7 3f  -16 

8.088b -23 

7 

2 . / 12b  -y6 

/ . 32ilf  -08 

0 . P 0 0 b - 0 1 

6 . “ 7 0 b - 1 6 

S.870b-|6 

1 .123b-22 

8 

2.”26b -06 

1 .G/SL-O/ 

0 . 00 oL -0 1 

3.020b -|6 

3.020b -16 

I .J7ib-22 

9 

3.008b  - Of) 

3.  39')b  -0  / 

0 . 0 0 0 1 - 0 1 

1 . 7SSt-i6 

1 . /SSL  - 1 6 

1 .9nHb-22 

1 0 

3 . 00  1 b -06 

7 . 32  1 b - 0 7 

0 . 0 (1  0 b "01 

1 . 0 3 7b  - 1 6 

1 . 0 3 71-16 

2. 06  lb  -22 

1 1 

0.20  Ob -06 

1 .S/SF-06 

0 . 000b -01 

0 . Oo6L- 1 7 

0 . 0/i6b  - 1 7 

2.9291-22 

I 2 

6 . 0 6 .1  b - 0 6 

3.  3'-'Sb  -06 

0 . Oil  Ob  - 0 1 

1 .006b  - 1 / 

1 .0061-1  7 

3.180F-22 

1 s 

0 . OPOb -06 

7 . 32ob -06 

0 . 0 0 0 b - 11  1 

8 .26 1 b - I 9 

'b.26lt-l9 

3.220E-22 

* * 

******* 

***** 

******* 

***** 

* * * 

* • f 

xj-  cur  inN  rfvnM 

riFV 
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Sample  Problem  6 


This  problem  is  also  similar  to  sample  problem  4 with 
the  exception  that  the  prompt  gamma  ray  source  is  a fission 
delta  function  source.  The  18  source  values  are  entered  via 
the  *G-SVAL  command. 


92 


TDATR  SAMPLE  PROBLEM  6 - INPUT 


‘TiTue 

•g-sval 

6. 

• EXC 


FIPStON  source  prompt  gamma  pays 

b.R7?-i> 

^ 3.R35-C’  T.F'So-p  233-2  2,116-?  /, ‘185-3 


1^^«1  .1'40|7 
4.23-3  6.79-u 


) no«2 
I .se-4 


.10729 
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TDATR  SAMPLE  PROBLEM  6 - OUTPUT 


f 

I 


TOATW  PROHLfM  NUMOtR  i KISSIUN  SOURCf  PROMPT  GAMMA  RAyS 

gamma  snoRCE  nurmA).  r/A r IMN=  i.oope^op  gammas/kt 

TlElO:  l.iliKlT+nO  KT  TMIAI..  ||IIPUT=  t.0n0f+0(l  gammas 

rtMf(SEr)  / T’'itijr.y  i ^^^Rl,y  ^ PNi-RGy  s entrgy  u energy  s energy  b 

O,n(l0Of-Ol  -Or’  1 . 5'l''E  -ll<?  'E.tAU^-nR  6.>W,’|-iTR  ><.h/«r-n?  t.7h«E-01 


I I ME  (SET)  / ENERGY  / energy  « E,N(R(.r  •)  ENERGY  Hi  ENERGY  II  ENERGY  \? 


0 . AonoE.-r)  1 i.'TiiPE-ni  l.noaE-oi  l.n;.n-oi  g-ihie-o.^  -^.r^Sf-o?  s./Sge-op 

<l»**1itt*lt*tltlttkttt********ll*‘:**k*** 

tlMllSEri  / ENEE’GY  lA  ENIRGY  pi  EhlRGY  IS  ENtR(;Y  1 (,  ENlRCY  17  ENERGY  18 
p.noooE-ni  r.iige-ar  3.R3nr_-i'3  g.yroe-p'J  i.G8oE-oa 

• *********«4-**««*»-*4****«*«**«***** 

T[)ArW  titlMMFR  s FISSIIVj  GAMMA  PAYS 

*«***«*««**«**««*****♦«** 


♦ * * * 


* * * * 


('.ROUND 

1 E VM 

n , (1  'in K M , 

n . iT(,  Ai.M/r ’!♦  * ,■> , 

•A  . n (1  0 N E T , 

n.onoMtlES 

*M.1R  1 1 . 

PA|jM 

A , M '1  n K M « 

MM  . R 3RI,M/(;  M‘ 

• GPSk  E T f 

n . '1  '<  7 M 1 1 E s 

SI  ANI 

rat  GE 

n _ ii  I'K  M , 

MG.'i  3Rr,M/C’'*  »R, 

M.sPS"T  E. 

'T.AR7MILE  S 

T AE  GE  1 

Al  T . 

M . 

1 .ItOoK’-, 

1 1 G , / (1 1 »,m/(;m * ♦? , 

3.r>H|AE  T, 

O.GRlMlLES 

siu)R(',r 

al  t . 

1 . I'AOIP  " , 

1 1 G . / pj  1 G.  M / r ki  ♦ ♦ p , 

1.?MIkPT, 

0 .GR  1 M I|.E  s 

♦ SI,  A N I 

AN'I.I  E 

A 

11  _ piroThE  1 

■ Ri.t  s ( r'is=  1 , oppoo  1 

*G  M 

olaTEh  erom  oiper  gop  iRD  I n a Te.-S 

ft  * * * 

♦ < > 

ft  ft 

ft  ft  ft  ft  ft 

***ftftftftft* 

ft  ft  « ft  * 

♦ ♦*♦♦♦ 

Bbljv 


V^i  i 
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TDATR  SAMPLE  PROBLEM  6 - OUTPUT  (Cont'd) 


pt-'nMPi  r.iiNfitML  nnsK  (wao/sfc) 


AUSOI  llr^  Ll.'CAL 


1 IMF 

r 1 MF 

"NFI'I  L IDEO 

SCAllEPFO 

TtiTaL 

cumulative 

1 

-00 

M.OlM'f'  0 

7 . ?9  SF  - ! 0 

0 . (lOOE- 

0 1 

/.2P5E-1'1 

7.2S5E-a5 

? . o7oi -no 

1 .s//F-np 

0 . OOOF -0  I 

0 . 1>0  0E  “ 

IS 

0 , so  oF  - I 5 

H.OSlE-25 

>, 

^.o7?E-00 

5.  x9/F-n9 

0 . Oflot  ”0  1 

S.050E- 

IS 

S,o  IbE-iS 

'>.US2t-25 

u 

(’.070E.  - flO 

7.5?i)E-oP 

0 , 0 n n 1 - 0 1 

5.601F- 

IS 

3.buiE-iS 

1 . iaoE-22 

5 

2 . OM'JE  - OO 

1 .S7SE-00 

0 . OOoE-O 1 

2.5PIF- 

|S 

2.321E-1S 

1 .007t-22 

4> 

7o<?E  -OO 

5.  JP'jF  -0« 

n.oooE'Oi 

1 .20F,f- 

IS 

1 .2obF-i5 

1 .72?.E-22 

f 

7U^E-00 

7.5cIiE-0O 

0 . nooF -0 1 

0,  /UI|^  •’ 

lO 

0,  7oyE-jO 

2.0P'1E-22 

fl 

a."?oF-oo 

t . S7SF -0 / 

o.nor.E-oi 

5. 5S7E- 

l^* 

5.5S7F-lb 

2.«7oE-22 

n 

5.00''E-nO 

i.  l9SE-n7 

n.nooE ”0 1 

t . 

lO 

1 . 75nE-iO 

2.‘»OOE-22 

10 

5.O0lf “OO 

7.5^oF-o7 

o.oooE-oi 

I.iaHE- 

1 0 

1 . 12«E-i6 

5.5nE-22 

1 1 

o.?'i'iF-no 

I .S7SF-PO 

n.nnoF-ni 

7. 1S2F- 

17 

7. 1S2E-17 

9,5nt-22 

l^ 

0 . 0 0 '1 E - n 0 

l'7SF-oO 

0 . oooE “0  1 

?.0OSF- 

1^ 

p.boSE-i 7 

a .‘7')oE-22 

15 

9.90<JF-ob 

'OO 

O.OOOE-Ol 

2. 50OE- 

1« 

2.504F-ie 

5. 12nE-22 

* 

* * « « ^ * Jt 

***** 

» * t * * * 

* * * * * 

* * 

***** 

A * * 

PkiIMPT  i;amM4  AfW  DUSK  (PAD/StC) 


AHSHLIITF 

time 

Lur  Ai 

1 1 “'t 

MUCFtELnif  f) 

scat  ifpfu 

TIUaL 

CIIMULAl  1 VE 

t 

2 . bhPF -Ob 

S . oi)H[  -to 

0 . PSb^ - 1 a 

0 . OOOF -n 1 

o.9SOF-ia 

0.95hF-25 

? 

2.m/I'F  - Ob 

1 .S7  /(  -OP 

n . oooF  -0 1 

0.220F -IS 

b.220E-1 S 

7 .0  7'iF  -2  5 

^ , r.  /Pi  -Ob 

Vi7I  -(,9 

0 .oonF -0 1 

S.  5/PF--1S 

S. 579E-1S 

P.Ol  |F -p5 

Q 

2 . b /OF  -Ob 

/.  52, IF  -oP 

0 . 0 0 n F - 0 1 

5 ,S{)  5F.  - 1 S 

5.Sn5E-lS 

1 .OOPt-p? 

s 

2.00'iE-ob 

1 .S  / tF  -oH 

0 .oooF -0 1 

2.2SOE-1S 

2.2SSF-1S 

1 .5'F0F-22 

h 

2 . 702> -Ob 

5 . 5PSF  -ciM 

0 . 0 0 0 F - 0 1 

1 . 2 a 0 F - 1 5 

1 .2anF-is 

1 .6S7F -pp 

1 

? . 70  2E-0b 

7 . ^20F ”0^ 

0 . 0 OoF  -0  1 

b.SPiE-tb 

h.S9|F-iF, 

P.OIOF-PP 

2. «2ol -no 

1 .S  /SF  -0  / 

0 . 0 n ii  F.  - 0 1 

5.  1 77E- 1 b 

5. 1 7 /F -lb 

P.  5/0F'P2 

1 

1 . 0 0 p t -no 

5.  <9')F  -0  / 

0 . 0 0 0 F - 0 t 

1 . a 1 «t  - 1 0 

1 . a 1 OF  - 1 b 

P.72'7F-pp 

\ 0 

5 . (j  0 1 F - 0 0 

7 . 52  0*  - 0 / 

o.oooE-ni 

s . 9 0 j E - 1 / 

S • ‘7  0 a F - 1 7 

5.0aSF -pp 

1 1 

a . 2 0 <i  f -Ob 

1 . S 7SF  -f|b 

0 . 0 0 0 r - 0 1 

2. 1 OPE- 1 7 

2. 1 opF-t  7 

5.PF'/F-P2 

0 , ohOF-OO 

5.  59SF -Ob 

0 . oooF -01 

S.aFFSE-|H 

S.'ifFSF-i« 

5.apaF -pp 

1 ^ 

9.9HPE-00 

7 . 520F  -Ob 

0 . 0 0 0 F - 0 1 

S.Fi.a  sE-t9 

i.««5F-l‘7 

5.  laaF-pp 

* 

*♦»»»*♦ 

***** 

****** 

****** 

****** 

* * * 

* • ( Xt  C'l  I IfiM  CIIXPI  M F'' 
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Sample  Problem  7 


The  geometry  and  response  configurations  for  this 
problem  are  the  same  as  for  sample  problem  4.  The  source  is 
a time  dependent  Gaussian  in  the  8-10  MeV  source  group.  The 
parameters  for  the  Gaussian  are:  the  mean  is  at  10  shakes 
the  amplitude  is  1 MeV/second  and  the  full-width-at-half- 
maximum  (FWHM)  is  two  shakes.  The  source  time  bins  are 
given  by  the  *G-TIME  command  in  shakes  (the  default  unit). 


TDATR  SAMPLE  PROBLEM  7 - INPUT 


♦title  GAllSSTA^4  SOURCt  PROMPT  GAMMA  RAYS 
*G-Sl)liRrF(2) 

♦ OSVAL  17*0  1 

*G-TIME  OP^lfe^RlOlliatOtbaOSO 

*G-RAIjSS  10  1 2 

*FXC 


97 


TDATR  SAMPLE  PROBLEM  7 - OUTPUT 


rOATW  PRf’HLFM  NUMHEH 


r.AllSSIAN  SftUKCF  PROMPT  GAMMA  RAYS 


****** 

« * « * A 

* * * 4 * 

* « A A A 

***** 

***** 

* * * * 

GahmA  $ni(RCE 

NORmAi  t /a  t lilN  = 

1^-01 

GAmmAS/K  t 

Y JKI  l)=  1 

.oooF+on 

KI  TMTAL  (1IJTPIJI  = 

a.3(,IF-0l 

GAMMAS 

TIMMSKC)  / 

ENF  RGY 

1 FNFRf.Y 

^ fnfri;y 

i energy  u 

1 ENFWGY  F 

1 ENERGY  h 

n,no<YOT-oi 

0 . 0 0 n K - 0 1 

P .POOF -P 1 

O.OPOF-Ol 

O.POPE-Pl 

P.POPF -01 

P.OOOF-Ol 

1? , 1)  (1 1'  0 T - n H 

0 .oonF -0 1 

O.oopF -01 

0 .pool -0 1 

0 .PPPF -P  1 

n. POOF-PI 

0 .OOpF -0  1 

« . oooof-  -nH 

0 , 0 p n F - 0 1 

O.POPF-Pl 

0 . 0 P P F - 0 1 

0 , n n 0 F - P 1 

P.OOPF-OI 

P.OOOF-01 

M.nonnh-nA 

p .opoE-o 1 

0 . OPOt-p  t 

P.POPF -0 1 

0 .POPF-O 1 

O.OOOF-OI 

O.POOF-Ot 

H.npnfiT  -n» 

O.OOOF-01 

O.POPL-Pl 

P.POPF -01 

0 .POPE-O 1 

P.POPF -01 

P.OOOF-Ol 

R.oon('t:-0'< 

p . 0 P 0 F - 0 1 

O.PPPF-Pl 

0. POPE-PI 

O.OOPF-Ol 

0 , OOOE-P 1 

0 ,00OL-0 1 

1 .nonoF.-o  / 

o.oopF  -m 

P.OPPF-m 

0 . 0 0 0 F - p 1 

p .POOF-0  1 

O.OOPF-01 

O.OOPF -0 1 

1 . TonOF -P  / 

0 . n (.1 0 F - 0 1 

O.PPOE-ol 

P.POPF-ol 

P.POPF -0 1 

O.OOOF-OI 

O.OOPF-01 

1 .2ni)nF.-o  1 

p.nnoF -01 

O.POpF.-fll 

P.OPPF-OI 

0 .POOF -n 1 

P . OPPF -0 1 

P .POOt >0  1 

1 .RnnoF-O  1 

0 . 0 p n F - 0 1 

P.OOpF-Pl 

P . 0 0 0 F - 0 1 

P.POPF -0 1 

O.POOE-Pl 

0 .onoE-0 1 

l ,60PiiF - 0 7 

O.OOPF-OI 

P.POPF -PI 

P.OPPF-01 

P.fiOOE-0  1 

p,onoE-oi 

O.OOOE-Ol 

^.OPOPF -0  7 

0 . OOoF -0 1 

0 . OOPF -p  1 

P.PPOF-01 

P.OOOE-Pl 

P .OOPF -0 1 

O.OOOE-01 

3.0000E-0  f 

o.nooF-ot 

P.POPl -01 

O.OOpF-oi 

P .POOE-O 1 

O.OPOF-01 

P.OOOF-Ol 

****** 

* * * * * 

* * * * * 

* * * * * 

* * » • • 

***** 

* * * * 

TiMKSfcCI  / 

fmf  wgv 

7 fmERGY 

R FNFRGY 

9 ENERGY  m 

• ENF  RGY  11 

INFRGY  1? 

O.OOPPt-01 

P.OPOF -01 

P.POPF -PI 

0.P0PF-01 

n.oppF-Pl 

P .POOF -0 1 

0,0POF-01 

P.onoijF-oB 

0 . PiipF  -p  ) 

P.popF -01 

0 .PPPF -0 1 

p.nooL-oi 

P . 0 P P F - 0 1 

P.OOOE-01 

a . onooF-OH 

o.ooor-(M 

P.POPl -01 

0,0001-01 

P.OPPF-Pl 

P.POPF -01 

P.OPOF -PI 

M.OPPPF-Pft 

0 , 0 (i  P F - U 1 

P . P P p F - p 1 

p.POPF-ni 

o.nnnF -oi 

0 .OppF -0 1 

P.PPOt-Ol 

M.OPOOF-PH 

O.OOPF -pi 

p . OOPF -0 1 

0 , OOOF -0 1 

P .OpoE -0 1 

0 .POOF “0 1 

P.POpE-01 

R.OPOOF  -n(il 

O.OOPF -01 

n . OPOF -p 1 

0 . 0 0 0 E - n 1 

n.oppF-oi 

0 . POPF -0 1 

P.OOOF-Ol 

, i.oopriF.-n7 

0 . PooF -0 1 

O.PPoF -01 

O.OOpF. -0  1 

P.OPOF-Pl 

P . P 0 0 F - fi  1 

0 .pool -P 1 

1 , 1 (iPi'F  -0  / 

o.poot-ni 

O.PPoF-Ot 

o.PooE-ni 

0 .POOE-p 1 

0 . nnoF -p 1 

P.OPOf-OI 

1 .PnopF -0  7 

0 . n 0 p F - 0 1 

0 . 0 P p F - 0 1 

0 . pppF -p 1 

0 ,npnt-o 1 

P . PPOF -p 1 

P.OOpt'Ol 

t .«nOOF -0  7 

o.nopt.-oi 

0 . 0 0 P F - 0 1 

0 . POpF -0 1 

p.nnpE-iii 

n . pppF -0 1 

O.OPpE-Ol 

T .MPPPF -0  7 

O.oopF -01 

P.POPl -01 

0 .OOpF -« 1 

o.pppE-ni 

P.PPOE-Pl 

0 .OpOt-0 1 

P,  OPPuF  -ft  7 

p . n 0 p F - p 1 

0 . popL-p 1 

0 .OOpF -0 1 

p ,onpF -0 1 

0 .OOPF -0 1 

P.OOPl-Ol 

^.OOOPE-pY 

P .OOPF -0 1 

P .POpF -P 1 

0 .OppF -p  1 

P .POpF -0 1 

P.POPF.-Ol 

P.OOPk-Ol 

****** 

* * * * ♦ 

* * * * * 

A A A A ♦ 

***** 

* * ♦ * * 

* * * • 

TDATR  SAMPLE  PROBLEM  7 - OUTPUT  (Cont'd) 


TIMf(StC)  ' 15  KNtPl.V  m I^NKRGV  is  KNt«GV  16  FNERGV  17  FNtt'Ct  18 


0. 0nnnp-ol  ti.ooni-oi  o.onoi-oi  o.oooF-m  o,ooot-oi  o.oooe”01  «.765t“52 

?.onnoF-o«  n.nont-OI  n.oont-oi  o.oonF-oi  o.O0oi-r>l  o.oooE-ot 

«.onnoK-o«  o.oonf-ol  o.ooni^-ni  o.ftooF-nt  o.oooF-oi  o.ooot-oi  i.6t/E-i2 

s. onooE-oa  n.onn^-ot  n.onoE-oi  o.onnE-oi  OtnooE-ol  o.onot-nt  i.69Sf-o6 

fl-ocnot-o^  n.oooE'i'l  o.finftE-ni  n.nooE“ni  o.oooE“nl  o.onoE-oi  6,9«of-o5 

9,ooooF-ofl  o.nonC-oi  o.ooot-oi  y.aooE-ni  o.^ooE-ni  o.oooF'Ol  S,sljfcF“02 

t. ooonE-07  o.oooE-oi  n.oooi-oi  o.nnoE-Pi  n.oooE-oi  o.noot-oi  i.niE-0| 

l.ir,noE-o7  o.onoE-ot  n.oooE-ol  o.onoE-nl  o.oooE-ni  o.onoE-oi  S.S5<>r-02 

l,2Pnof.-p7  o.oonE-oi  o.oooE-oi  o,oonE-(n  O.OOPE-Ot  o.nooE-oi  s.vaaE-oj 

1. <»0n0F-O^  n,OOnF-Ol  O.OOOE-Ot  O.OnnE-ol  O.noOF-OI  O,O00E-01  1.6<J5t-06 

l.^>oooE-o7  n.onni-01  o.oooF-oi  o.nnot-ol  p.onot-oi  o«oooF-ot  t.6t7E-i^ 

2.00no6-o7  0.00(lfc-01  O.OOOE-01  O.npnK-ol  0,ilOOf-Ol  0*nOOF-01  «.763F-12 

3.000i)E-o7  n.OOoF-Ot  O.OOOF-Ol  O.OOOE-Ol  o.onoE-ol  o.ooot-oi  O.OOOF-fl 
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